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Abstract 

We review the study of inhomogeneous perturbations about a homogeneous and 
isotropic background cosmology. We adopt a coordinate based approach, but give 
geometrical interpretations of metric perturbations in terms of the expansion, shear 
and curvature of constant-time hypersurfaces and the orthogonal timelike vector 
field. We give the gauge transformation rules for metric and matter variables at 
first and second order. We show how gauge invariant variables are constructed by 
identifying geometric or matter variables in physically-defined coordinate systems, 
and give the relations between many commonly used gauge- invariant variables. In 
particular we show how the Einstein equations or energy-momentum conservation 
can be used to obtain simple evolution equations at linear order, and discuss exten- 
sions to non-linear order. We present evolution equations for systems with multiple 
interacting fluids and scalar fields, identifying adiabatic and entropy perturbations. 
As an application we consider the origin of primordial curvature and isocurvature 
perturbations from field perturbations during inflation in the very early universe. 
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1 Introduction 



The standard model of hot big bang cosmology is based upon the spatially 
homogeneous and isotropic Friedmann-Robertson- Walker (FRW) model. This 
successfully describes the average expansion of the universe on large scales 
according to Einstein's theory of general relativity, and the evolution from a 
hot, dense initial state dominated by radiation to the cool, low density state 
dominated by non-relativistic matter and, apparently, vacuum energy at the 
present day. This standard model is described by just a handful of numbers 
specifying the expansion rate, the temperature of the present microwave back- 
ground radiation, the density of visible matter, dark matter and dark vacuum 
energy. 

But an homogeneous model cannot describe the complexity of the actual distri- 
bution of matter and energy in our observed universe where stars and galaxies 
form clusters and superclusters of galaxies across a wide range of scales. For 
this we need to be able to describe spatial inhomogeneity and anisotropy. But 
there are few exact solutions of general relativity that incorporate spatially in- 
homogeneous and anisotropic matter and hence geometry. Therefore we use a 
perturbative approach starting from the spatially homogeneous and isotropic 
FRW model as a background solution with simple properties within which we 
can study the increasing complexity of inhomogeneous perturbations order by 
order. 

The introduction of a homogeneous background model to describe the inho- 
mogeneous universe leads to an ambiguity in the choice of coordinates. In the 
FRW model the homogeneous three-dimensional hypcrsurfaces provide a nat- 
ural time slicing of four-dimensional spacetime. For instance, hypersurfaces of 
a uniform density coincide with hypersurfaces with uniform spatial curvature. 
However in the real, inhomogeneous universe spatial hypersurfaces of uniform 
density do not in general have uniform spatial curvature, and hypersurfaces of 
uniform curvature do not have uniform density. In general relativity there is, 
a priori, no preferred choice of coordinates. Choosing a set of coordinates in 
the inhomogeneous universe which will then be described by an FRW model 
plus perturbations amounts to assigning a mapping between spacetime points 
in the inhomogeneous universe and the homogeneous background model. The 
freedom in this choice is the gauge freedom, or gauge problem, in general rel- 
ativistic perturbation theory. It can lead to apparently different descriptions 
of the same physical solution simply due to the choice of coordinates. This 
freedom can be a powerful tool as it allows one to work in terms of variables 
best suited to the problem in hand, however it is important to work with 
gauge-invariant definitions of the perturbations variables if the results are to 
be easily assimilated and cross-referenced to the work of others. 
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In this review we will focus on how one can construct a variety of gauge- 
invariant variables to deal with perturbations in different cosmological models 
at first order and beyond. We will emphasise the geometrical meaning of metric 
and matter perturbations and their gauge-invariant definitions. 

Most work to date has been done only to linear order where the perturbations 
obey linear field equations. Even then we must consistently solve the linear 
evolution equations, subject to the constraint equations of general relativity 
which relate the matter variables to the geometry. Beyond first order, the non- 
linearity of Einstein's equations becomes evident, making progress much more 
difficult. But in limiting cases, notably the large scale limit, it is possible to 
extend some of the simple results of linear theory to higher orders. 

Previous reviews on the topic of linear perturbation theory in cosmology in- 
clude, for example, Refs. [66,123,43,21,141,61], and the relevant chapters in 
Refs. [177,131,75,84]. For second-order perturbation theory see Refs. [129,15,126] 
For reviews on inflation in particular see, e.g., Refs. [88,70,86,93,84,77,17,174] 
and for cosmic microwave background anisotropies see, e.g., Refs. [67,82,42,51] 
Finally, for reviews on perturbation theory in the context of the higher- 
dimensional brane-world scenario see, e.g., Refs. [68,24,135,100,23] and for 
perturbation theory in the context of the low-energy string effective action 
see, e.g., Refs. [85,34]. 

For simplicity we work with a flat background spatial metric which is com- 
patible with current observations. For generalisation to spatially hyperbolic 
or spherical FRW models sec other reviews, e.g., Rcf.[66]. We use a prime to 
denote derivatives with respect to conformal time, and we use a comma to 
denote partial derivatives with respect to comoving spatial coordinates, i.e., 

T = — 

For further definitions and notation see Appendix A. 
2 Perturbations in cosmology 

Throughout this review we will assume that our observable Universe can ap- 
proximately be described by a homogeneous and isotropic Friedmann- Robertson- 
Walker (FRW) spacetime. Thus we assume that the physical quantities can 
usefully be decomposed into a homogeneous background, where quantities 
depend solely on cosmic time, and inhomogeneous perturbations. The pertur- 
bations thus "live" on the background spacetime and it is this background 
spacetime which is used to split four- dimensional spacetime into spatial three- 
hypersurfaces, using a (3+1) decomposition. In addition we work with a spa- 



(1.1) 
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tially flat FRW background, though results can be readily extended to non- 
zero spatial curvature. 

We start this section by defining arbitrary perturbations of tensorial quantities 
and then proceed by decomposing vectors and tensors into "time" and "space" 
parts on the spatial hypersurfaces. 

2.1 Defining perturbations 

Any tensorial quantities can then be split into a homogeneous background and 
an inhomogeneous perturbation 

T{ri,x')=To{ri) + 5T{ri,x'). (2.1) 

The background part is a time-dependent quantity, Tq = To(?7), whereas 
the perturbations depend on time and space coordinates x'^ — [t],x'^]. The 
perturbation can be further expanded as a power series, 

oo n 
n=l ^■ 

where the subscript n denotes the order of the perturbations, and we explicitly 
include here the small parameter e. In linear perturbation theory, for example, 
we only consider first-order terms, e^, and can neglect terms resulting from 
the product of two perturbations, which would necessarily be of order or 
higher, which considerably simplifies the resulting equations. In the following 
sections we shall omit the small parameter e whenever possible, as is usually 
done to avoid the equations getting too cluttered. 

2.2 Decomposing tensorial quantities 

It is convenient to slice the spacetime manifold into a one-parameter family 
of spatial hypersurfaces of constant time, which is the standard (3-1-1) split 
of spacetime. This foliation was introduced by Darmois already in 1927 (see 
Ref. [53]) and popularised by Arnowitt, Descr and Misncr [6] (for conditions 
on the existence of the foliation see e.g. Ref. [172]). We refer to the foliation 
of spacetime by spatial hypersurfaces of given conformal time as the time- 
slicing, and the identification of spatial coordinates on each hypersurface as 
the threading. 

Note that the (3-1-1) split of spacetime precedes the decomposition of 3- 
dimensional quantities into scalars and vectors, or scalar, vectors, and tensors 
below. 
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2.2.1 Vectors 



We can split any 4 vector into a temporal and spatial part, 



W^W' . (2.3) 



Note that is a scalar on spatial hypersurfaces. The spatial part W can then 
be further decomposed into a further scalar part U and a vector part U^^^, 

W = S'^Uj+Ui^, (2.4) 

where dU^^/dx'' = 0. The denominations scalar and vector part go back 
to Bardeen [7] and are due to the transformation behaviour of U and Ul^^ 
on spatial hypersurfaces. The decomposition of a vector field into potential 
(or curl-free) part and a divergence- free part in Euclidean space is known as 
Helmholtz's theorem. The curl-free and divergence-free parts are also called 
longitudinal and solenoidal parts, respectively. 

In our isotropic (FRW) background, there can be no spatial vector part at 
zeroth order (as this would correspond to a preferred direction), but there can 
be a non-zero temporal part, 

WfoT^O, Wjo) = 0. (2.5) 

A non-trivial vector, with a non-zero spatial part, can appear only at first 
order. We give examples of vector fields in a FRW background including per- 
turbations in Section 3.1, where we discuss the unit vector field normal to 
constant-time hypersurfaces and the fluid 4-velocity. 

Note that "divergence-free", etc, is defined with respect to the fiat-space met- 
ric, rather than using covariant derivatives, since perturbations are defined 
with respect to the spatially flat background. 



2.2.2 Tensors 

As for vector flelds, we can decompose a rank-2 tensor into a time part and 
spatial part, but now have also mixed time-space parts. 

Consider the metric tensor gi^i^ which we require to be symmetric, 

Qfiu = 9v,x ■ (2.6) 

The metric tensor has therefore only 10 independent components in 4 dimen- 
sions. We first split the metric tensor into a background and a perturbed part, 
using Eq. (2.1). It then turns out to be useful to split the metric perturbation 
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into different parts labelled scalar, vector or tensor according to their trans- 
formation properties on spatial hyp ersurf aces [7,162], which are themselves 
expanded into first and higher order parts using Eq. (2.2). 

Our background spacetime is described by a spatially flat FRW background 
metric 

ds^ = -drj^ + 5ijdx'dx^ , (2.7) 

where rj is conformal time and a = a{ri) is the scale factor. The cosmic time, 
measured by observers at fixed comoving spatial coordinates, x\ is given by 
t = /«('?) df]- 

The perturbed part of the metric tensor can be written as 



Sgoo^-2a^cf>, (2.8) 
Sgoi = a'^Bi , (2.9) 
Sgij^2a^Cij. (2.10) 

The — i and the i — j-components of the metric tensor can be further de- 
composed into scalar, vector and tensor parts 



Bi^B,i-Si, (2.11) 
Q. = -^/, + E,ij + F(,,,) + . (2.12) 

where 0, B, ip and E are scalar metric perturbations. Si and Fi are vector 
metric perturbations, and hij is a tensor metric perturbation, which we will 
now define. 

S'ca/ar perturbations can always be constructed from a scalar, or its derivatives, 
and background quantities. Any 3- vector, such as B i, constructed from a scalar 
is necessarily curl-free, i.e., B^[ijj = 0. 

Vector perturbations are divergence- free. For instance one can distinguish an 
intrinsically vector part of the metric perturbation 6goi, which we denote by 
—Si, which gives a non- vanishing contribution to Sg^^ij]. Similarly we define 
the vector contribution to 5gij constructed from the (symmetric) derivative of 
a vector F(i j) . 

Finally there is a tensor contribution to 6gij = a^hij which is both transverse, 
hij^ = (i.e., divergence-free), and trace- free {h\ = 0) which therefore cannot 
be constructed from inhomogeneous scalar or vector perturbations. 

Note that when raising and lowering spatial indices of vector and tensor per- 
turbations we use the comoving background spatial metric, 5ij, so that, for 
instance, hi = 5^^hik. 
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The reason for splitting the metric perturbation into these three types is that 
the governing equations decouple at linear order, and hence we can solve each 
perturbation type separately. At higher order this is no longer true: we find for 
example at second order that although the "true" second order perturbations, 
^g2ixvi still decouple, but their governing equations have source terms quadratic 
in the first order variables, 5gi^v. mixing the different types [127]. Indeed at all 
higher orders, n > 1, the different types of perturbations of order n decouple, 
but are sourced by terms comprising perturbations of lower order. 

We have introduced four scalar functions, two spatial vector valued functions 
with three components each, and a symmetric spatial tensor with six compo- 
nents. But these functions are subject to several constraints: hij is transverse 
and traceless, which contributes four constraints, Fi and Si are divergence-free, 
one constraint each. We are therefore finally left 10 new degrees of freedom, 
the same number as the independent components of the perturbed metric. 

The choice of variables to describe the perturbed metric is not unique. Already 
at first order there are different conventions for the split of the spatial part of 
the metric. We follow the notation of Mukhanov ct al [123] so that the metric 
perturbation, ip, can be identified directly with the intrinsic scalar curvature 
of spatial hypersurfaces at first order, see later. Sometimes it is useful to work 
in terms of the trace of the perturbed spatial metric 

C = Q = + V^^; (2.13) 

At first order this coincides with the perturbation of the determinant of the 
spatial metric. Including terms up to second order we have 

det + 2Cij) = 1 + 2C + 2 (C^ - dfi'^) 

= 1 - 6^ + 2V'^E 
+12^^^ - Si^V'^E + 2{V'^Ef - 2E,ijE/^ - 2FijF\^ 

~hij0 - 2E^ij0 - 2Fij0 . (2.14) 

where we have used the general result det(e'^) = e'^^'''^. There are further 
choices for the way the spatial metric is split into the different perturbation 
variables at second (and higher) order in the perturbations. 

Note that our metric perturbations in Eqs. (2.8-2.12) include all orders. If we 
write out the complete metric tensor, up to and including second-order in the 
perturbations we have 
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^00 = -a' (1 + 201 + 02) , 
Qij = [Sij + 2Ciij + C2ij] , 



(2.15) 



where the first and second order perturbations Bu and Cuj, and and C2ij, 
can be further split according to Eqs. (2.11) and (2.12). 

The contravariant metric tensor follows from the constraint (to the required 
order) , 

g^^ug'^^S^^ (2.16) 

which up to second order gives 



,00 . 



a 

-2 



201 



^2 + 



BikB^ 



g'^^a-' 



Bi + -B2^ — 20iSJ — 2SifcCi*^' 



(2.17) 



3 Geometry of spatial hypersurfaces 



3.1 Timelike vector fields 



The perturbed metric given in Section 2.2.2 implicitly defines a unit time-like 
vector field orthogonal to constant-r; hypersurfaces, 

oc ^ , (3.1) 

subject to the constraint 

n% = -1 . (3.2) 

In the FRW background this coincides with the 4- velocity of matter and the 
expansion of the velocity field 9 = 3H, where H is the Hubble expansion rate. 
We define the conformal Hubble parameter 

n = aH. (3.3) 



In this section we calculate corresponding geometrical quantities for n"^, and 
thus the space-time, defined by the perturbed metric tensor. Note however 
that in the perturbed spacetime the vector field n'^, need no longer coincide 
with the 4- velocity of matter fields at first order and beyond. 
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Up to and including to second order, the covariant vector field is 



= -a 



(3.4) 



and the contravariant vector field is 



n 



i[l-0,-i^ + 50?-iiS:.i3f 

ir i_\ « 



a 



(3.5) 



Observers moving along the hypersurface-orthogonal vector field, n^, have a 
vanishing 3-vclocity with respect to the spatial coordinates when the shift 
vector 5* is zero. We will refer to these as orthogonal coordinate systems. 



3.2 Geometrical quantities 



The covariant derivative of any time-like unit vector field can be decom- 
posed uniquely as follows [172]: 



where the spatial projection tensor V^^i orthogonal to n'', is given by 



(3.6) 



(3.7) 



The overall expansion rate is given by 

the (trace-free and symmetric) shear is 

1 



the (antisymmetric) vorticity is 



(3. 



(3.9) 



and the acceleration is 



(3.10) 



(3.11) 

On spatial hypersurfaces the expansion, shear, vorticity and acceleration co- 
incide with their Newtonian counterparts in fiuid dynamics [59,160]. In this 
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subsection wc focus on the unit normal vector field n'^, but the expansion, 
shear, vorticity and acceleration defined in this way can readily be applied 
to any other 4-vector field, such as the 4-velocity u'^. One can easily verify 
that the vorticity (3.10) is automatically zero for the hypersurface orthogonal 
vector field, njj, defined in Eq. (3.4). Note however that the perturbed fiuid 
velocity can have vorticity and this is described by the vector (divergence-free) 
part of the fiuid 3- velocity which we will define in section 4. 

The projection tensor V^i, is the induced 3-metric on the spatial hypersurfaces, 
and the Lie derivative of Vfj,i, along the vector field is the extrinsic curvature 
of the hypersurface embedded in the higher-dimensional spacetime [172,40]. 
The extrinsic curvature of the spatial hypersurfaces defined by 71^ is thus given 



At first order we can easily identify the metric perturbations with geometrical 
perturbations of the spatial hypersurfaces or the associated vector field, n^. 

The intrinsic curvature of spatial hypersurfaces up to first order is given by 



by 




(3.12) 



(3.13) 



a 



The scalar part of the shear (3.9) up to first order is given by 




(3.14) 



where we define the shear potential 



(3.15) 



The vector and tensor parts are given by, respectively. 





(3.16) 



The acceleration up to first order is 



(3.17) 



The expansion rate up to first order is given by 



a I 3 



(3.18) 



9 



The intrinsic spatial curvature, shear and acceleration of are given up to 
second order in Appendix C. 

The overall expansion, up to second order is given by 



(Jj (Jj 

-l^ (02 - 30,^) + I - B,^) + 0, - C,,'^') 



2 a 

3 a' 
'2^ 



(3.19) 



Focusing for the moment only on scalar perturbations (neglecting first order 
vectors and tensors) on large scales (neglecting spatial derivatives) the per- 
turbed part of the expansion simplifies to 



I 





fa' 


a 


^ a 


3 


'a' 


a 


a 



a 9 

3-0 
a 



(3.20) 
(3.21) 



Note that the expansion rate, 9 in Eq. (3.19), is defined with respect to proper 
time (comoving with n^). The expansion rate with respect to conformal time 
on large scales is 



e, 



coord 



(3.22) 



We see that on large scales (and considering only scalar perturbations) spa- 
tially flat hypersurfaces (^ = 0) are also uniform coordinate expansion hyper- 
surfaces (on which the perturbed expansion vanishes). 

In the following we will often refer to the perturbed logarithmic expansion, or 
"number of e-foldings" 



N 



1 f 1 r 1 

-J edr^ - J OcooTddr] = In a - V'l - -ip2 - ipl ■ (3.23) 



This quantity, and its perturbation 6N becomes a particularly useful quantity 
to describe the primordial scalar perturbation beyond linear order, as we will 
discuss in Section 11. 
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4 Energy-momentum tensor for fluids 



We are interested in how the spacetime geometry, described by the metric 
tensor, is affected by the perturbed matter content, described by the energy- 
momentum tensor. 

The four-velocity of matter, u^, is defined by 



(4.1) 



where r is the proper time comoving with the fluid, subject to the constraint 



(4.2) 



The components of the 4- velocity up to second order are then given by 



uq — —a 



1 + 01 + ^02 - ^0l' + ^Vikv'l 



Ui — a 



vii + Bu + - (v2i + B2i) - (t)iBu + 2Cukv'l 



u —a 



1 - 01 - -02 + -01 + -vikvt + VikBl 



(4.3) 



(4.4) 



The spatial part of the velocity can be split into a scalar (potential) part and 
a vector (solenoidal) part, order by order, following Eq. (2.4) as 

^ 5^^v, + , (4.5) 

for each order n where we refer to V(^n^ as the scalar velocity potential, and to 
'^ln)vec ^ vector part. 

Note that is the 3-velocity of matter defined with respect to the spatial 

coordinates, x^, and so is not the velocity with respect to the hypersurface- 
orthogonal vector field n^, except in orthogonal coordinate systems for which 
5* = 0. In comoving orthogonal coordinates, which we will discuss later, = 
and B' = 0. 
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4-1 Single fluid 



The energy-momentum tensor of a fluid with density p, isotropic pressure P 
and 4-velocity u^, defined above in Eq. (4.3), is defined as [66,160,47,177] 

T^,^ = {p+ P) u^u^ + Pg^y + TT^^ . (4.6) 

The anisotropic stress tensor tt^,^ is spht into first and second order parts in 
the usual way, 

and is subject to the constraints 

7r^,M" = 0, 7r'; = 0. (4.8) 

The anisotropic stress vanishes for a perfect fluid or minimally coupled scalar 
fields, but may be non-zero in the presence of free-streaming neutrinos or a 
non-minimally coupled scalar field. 

We follow Kodama and Sasaki [66] in defining the proper energy density as 
the eigenvalue of the energy-momentum tensor, and the four velocity as 
the corresponding eigenvector 



T^X = -pu" . (4.9) 



The anisotropic stress tensor decomposes into a trace-free scalar part, 11, a 
vector part, Ilj, and a tensor part, Uij, at each order according to 



2 

TTij — a 



(4.10) 

We get for the components of the stress energy tensor in the background 

n = -po, ^0=0, T'.^S'^Po, (4.11) 
at first order, 

W6Tl = -6p,, (4.12) 

(i)5T° = (po + Po)(t'i. + 5iO , (4.13) 

^'^ST.^5Pi5'j + a-\l)j, (4.14) 

and at second order 
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(2)5T° = (po + Po) 



(4.15) 



V2i + B2i + 4CukVi - 24)1 {vu + 2Bii) 

2 



+2 (5pi + 5Pi) (^;h + Bu) + - (b1 + vl) TTHfe , (4.16) 

a 

^''^ST. = 5P2 5^ + a-\' - \c'^^i3k + 2 (po + ^o) v\ {v,, + B,,) . (4.17) 



Note, that for compactness of presentation we have not split perturbations 
into their constituent scalar, vector and tensor parts in the above expressions. 
The decompositions are given for Vi and Tiij in Eqs. (4.5) and (4.10), and given 
for Bi and Cij in Section 2.2.2, in Eqs. (2.11) and (2.12). 

Note that contracting the i—j part of the energy-momentum tensor, Eq. (4.17), 
including the constraints for the anisotropic stress, Eq. (4.8), guarantees that 
the anisotropic stress cancels on the trace. 



W5r^fc = 35Pi, (4.18) 
^''^5T\ = 35P2 + 2 (po + ^o) v1 {v^k + Sife) . (4.19) 

This cancellation is true at all orders. 

Coordinate transformations affect the spUt between spatial and temporal com- 
ponents of the matter fields and so quantities like the density, pressure and 
3- velocity are gauge-dependent, as described in Section 6. Density and pressure 
are 4-scalar quantities which transform as given in Eq. (6.18) in the following 
section, but the 4-velocity is a 4-vector which transforms described in Sec- 
tion 6.3. The anisotropic stress is gauge-invariant at first order, but becomes 
gauge-dependent at second order. 



4.. 2 Multiple fluids 



In the multiple fluid case the total energy-momentum tensor is the sum of the 
energy-momentum tensors of the individual fluids 

T^'' ^Y.Tl'a)- (4-20) 

a 

For each fluid we deflne the local energy-momentum transfer 4-vector Q'^^) 
through the relation 

VMT(r) = QU , (4.21) 

where the energy- momentum tensor, T^'^'^, is locally conserved only for non- 
interacting fluids, for which Q"^^-^ — 0. Equations (8.2) and (4.21) imply the 
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constraint 



a 



(4.22) 



Following Refs. [66,111] we split the energy-momentum transfer 4- vector using 
the total fluid velocity as 

QU = Qi^)^' + /&) > (4-23) 

where Q[a) is the energy transfer rate and /^'^^ the momentum transfer rate, 
subject to the constraint 

n,fU = ' (4-24) 
From Eq. (4.24) we find for the temporal component of the momentum transfer 
rate vector fj^^^ 

= , = 2/^(„) (v,k + B,,) , (4.25) 

We then find for the temporal components of the energy transfer 4-vector to 
be at zeroth, first and second order, respectively 



la 



QU = ^ [^Q2a + Qoa (30l' - </>2) - 20i5Q 

/2 

+ (^^Ifc + B^k) I -AV) + QoaV^^ 



(4.26) 
(4.27) 

(4.28) 



where Qoa, ^Qia-, and 5Q2a are the energy transfer to the a-fiuid in the back- 
ground, at first and at second order, respectively. 

For the spatial components of the energy transfer 4-vector, the momentum 
part, we get at first and second order, respectively 



Q\o) — -QootV\ + —^fl(a) > 

-/2(a) + ^QlaV\ + Q,a (^^2 + 2<^li?l - ^C^A) 



1 

2a 



(4.29) 
(4.30) 



where /^^^ and /gj^^^ are the spatial parts of the momentum transfer rates at 
first and second order. 

Note that the homogeneous and isotropic FRW background excludes a zeroth 
order momentum transfer. 
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Using Eq. (2.4) the spatial momentum transfer vector of order n, /^q), can 
be further decomposed into a scalar and a vector part, according to 

fnia) ^ S''Ua),J + fn^a) ■ (4-31) 



5 Energy-momentum tensor for scaleir fields 

5.1 Single field 

A minimally coupled scalar field is specified by the Lagrangian density 

C = -lg'''<P,,<P,. - U{<p) , (5.1) 

where the scalar field kinetic energy is then non-negative for our choice of 
metric signature. 

The energy momentum tensor is defined as 

dC 

and we therefore get for a scalar field ip 

n = 9"'^,a^,. - K (U{^) + ^^'^Va^.a) . (5.3) 

Comparing to the energy-momentum tensor for a perfect fluids, Eq. (4.6), we 
can identify the non-linear 4-velocity, density and pressure of the scalar field 
[165] 

P = V>,A<^,« - • (5.4) 

Note that the anisotropic stress, n^u, is identically zero for minimally coupled 
scalar fields. In fact we can subdivide the energy-momentum tensor for a 
single field into that of a stiff kinetic fluid with p(<^) = P(<^) = — gf^'^(^^A¥',/t and 
a vacuum energy = —P(u) = U{ip), which exchange energy-momentum 
% = -Qiu) = {dU/dcp)V^ip. 

Sphtting the scalar field into a homogeneous background field and a pertur- 
bation, 

ip{ri,x') = ipo{ri) + 6(pi{rj,x') , (5.5) 
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and using the definitions above we find for the components of the energy 
momentum tensor of a perturbed scalar field at linear order without specifying 

a gauge yet 



7^0 = -2«~Vo' -Uo + a-^ip'o (01 ip'o - Vi) - U^Sipi , 

1 
2 



3 



-a Vo^ - Uq - U^^Sipi + a ^ip^ {dip[ - 4>i (p'o) 



3 ' 



(5.6) 
(5.7) 
(5.8) 



where U^^ = dU/dip and Uq = U{ipo). By comparing Eq. (5.8) with Eq. (4.14), 
we see that scalar fields do not support vector or tensor perturbations to first 
order. 



5.2 Multiple fields 



For N minimally coupled scalar fields the Lagrangian density is given by 



1 



(5.9) 



The energy-momentum tensor is 



T, 



N 

E 

K=l 



af3 



(5.10) 



where (pK is the Kth scalar field and U the scalar field potential and ipK,ij. = 

Analogous to the energy-momentum tensor for a single field, we can identify 
the non-hnear 4- velocity, density and pressure of each of the scalar fields [111] 



P{i) = P{i) = -g^>i,xfi,K ■ (5.11) 

The energy-momentum transfer to each fluid is = {dU/dipi)ipi^n. 

The total energy-momentum tensor (5.10) is the sum over N stiff fluids plus 
the vacuum energy 



N 



(5.12) 



1=1 
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The anisotropic stress, vr^jy, is identically zero for any number of minimally 
coupled scalar fields. The total energy-momentum tensor is only equivalent to 
that of a single scalar field in the special case where all the 4- velocities of the 
fields, are identical. This is true in the homogeneous FRW cosmology, 

but in general breaks down when inhomogeneous perturbations are considered. 

We split the scalar fields ipi into a background and perturbations up to and 
including second order according to Eqs. (2.1) and (2.2), 

Vi{.x^) = miv) + S^uix'^) + l5^2i{x'') . (5.13) 
The potential U = U {(pi) can be spht similarly according to 

U{ipi)^Uo + SUi + ^5U2, (5.14) 

where 

SUi = Y,UKSipiK, (5.15) 

K 

SU2 = ^UklS(PikS^Pil + J2^kS(P2K- (5.16) 

K,L K 

and we use the shorthand Uk = dU/d(pK- The energy-momentum tensor, 
Eq. (5.10), expanded up to second order in the perturbations for the metric 
tensor Eq. (2.15) is given in Appendix B.2. 



6 Gauge transformations 

We now review how tensorial quantities change under coordinate transforma- 
tions [7,162,121,26,108] (see Ref. [142,177,163] for earlier work on this subject). 
While the order of the perturbation is indicated by a subscript, we also keep 
the small parameter e in the following equations whenever appropriate. 

A problem which arises in cosmological perturbation theory is the presence 
of spurious coordinate artefacts or gauge modes in the calculations. Although 
the gauge modes had been dealt with on a "case by case" basis before, the 
gauge issue was resolved in a systematic way by Bardeen [7] . The gauge issue 
will arise in any approach to GR that splits quantities into a background and 
a perturbation. Although GR is covariant, i.e. manifestly coordinate choice 
independent, splitting variables into a background part and a perturbation is 
not a covariant procedure, and therefore introduces a coordinate or gauge de- 
pendence. By construction this only affects the perturbations; the background 
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quantities remain the same in the different coordinate systems. Here we as- 
sume the difference between the coordinate systems is small, of 0(e), however 
the gauge problem would persist also for finite transformations. Note that the 
"covariant approach" [46] also corresponds to a choice of gauge, namely the 
comoving one, which is made explicit by the inclusion of the velocity field 
[35,78]. 

In order to restore covariance as far as possible, we usually wish to eliminate 
the gauge degrees of freedom. We will show in Section 7 how, by construct- 
ing variables corresponding to perturbations in physically defined coordinate 
systems, the gauge dependencies can be made to cancel out (the quantities so 
constructed will not change under a gauge transformation). 



6.1 Active and passive approaches to gauge transformations 



There are two approaches to calculate how perturbations change under a small 
coordinate or gauge transformation. For the active view we study how pertur- 
bations change under a mapping, where the map directly induces the transfor- 
mation of the perturbed quantities. In the passive view the relation between 
two coordinate systems is specified, and we calculate how the perturbations 
are changed under this coordinate transformation. In the passive approach the 
transformation is taken at the same physical point, whereas in the active ap- 
proach the transformation of the perturbed quantities is evaluated at the same 
coordinate point. We will discuss both approaches briefly in the following, but 
shall use the active approach to calculate the transformation behaviour of the 
first and second order variables. For a mathematically more rigorous discussion 
see e.g. Ref. [108]. 



6.1.1 Active approach 

The starting point in the active approach is the exponential map, that allows 
us to immediately write down how a tensor T transforms up to second order, 
once the generator of the gauge transformation, has been specified. The 
exponential map is 

T = e-^sT, (6.1) 

where £^ denotes the Lie derivative with respect to The vector field gen- 
erating the transformation, is up to second order 

e ^ + le'^^ + 0{e') . (6.2) 
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The exponential map can be readily expanded 

exp(£^) = 1 + e£^, + ^6^4 + ^-e' + . . . (6.3) 

where we kept terms up to O(e^). Splitting the tensor T up to second order, as 
given in Eq. (2.2), and collecting terms of like order in e we find that tensorial 
quantities transform at zeroth, first and second order, respectively, as [121,26] 



To = To, (6.4) 
e(5Ti = e5Ti + ei:^,To, (6.5) 
e'5T2 = e' (8T2 + i^^To + 4,To + 2£^/Ti) . (6.6) 

Note that the background quantities arc not affected by the mapping. Wc will 
apply Eqs. (6.5) and (6.6) to scalars, vectors, and tensors after discussing the 
passive approach next. 

Applying the map (6.1) to the coordinate functions we get a relation for 
the coordinates of a point q in and a point p as 

x^'{q) = v x^ip) , (6.7) 

where we have used the fact that when acting on scalars £^ = ^^{d/dx^) and 
the partial derivatives are evaluated at p. The left-hand-side and the right- 
hand-side of Eq. (6.7) are evaluated at different points. Equation (6.7) can 
then be expanded up to second-order as 



x'^iq) = x'^ip) + eCnp) + \e' {iUP)ii(P) + ^^{p)) ■ (6-8) 

Note that we do not need Eq. (6.8) to calculate how perturbations change 
under a gauge transformation in the active approach, it simply tells us how 
the coordinates of the points p and q are related in this approach. 



6.1.2 Passive approach 

In the passive approach we specify the relation between two coordinate systems 
directly, and then calculate the change in the metric and matter variables when 
changing from one system to the other. As long as the two coordinate systems 
are related through a small perturbation, the functional form relating them is 
quite arbitrary. 

However, in order to make contact with the active approach, discussed above, 
we take Eq. (6.8) as our starting point. Note, that all quantities in the passive 
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approach are evaluated at the same physical point. Equation (6.8) can be 
rewritten to give a relation between the "old" (untilded) and the "new" (tilde) 
coordinate systems [26,108], 

x^{q) = x^{q) - + {it{q)A'^{q) - e2{q)) , (6.9) 

evaluated at the same physical point q. 

The passive point of view is very popular at first order, see e.g. the original 
paper by Bardeen [7] and the widely used reviews by Kodama and Sasaki [66], 
and Mukhanov, Feldman, and Brandenberger [123]. 

The starting point in the passive approach is to identify an invariant quantity, 
that allows us to relate quantities to be evaluated in the two coordinate sys- 
tems. We denote the two coordinate systems by and x^ system, and their 
relation is given by Eq. (6.9). We choose as an example the energy density, 
p, which as a four scalar will not change under a coordinate transformation. 
However, once it has been split into the background quantity and perturbation 
at different orders, these variables will change. 

The two coordinate systems are related by Eq. (6.9), which we can use to 
expand the right-hand-side of Eq. (6.11) in a Taylor expansion up to second 
order. To first order, the two coordinate systems are simply related, using the 
linear part of Eq. (6.8), by 

i^ = a;'^-^f. (6.10) 

We get the transformation behaviour of the density perturbation, 5p, from the 
requirement that the total density, p — po + Sp, has to be invariant under a 
change of coordinate system and therefore has to be the same in the x^^ and 
the x'' system, that is 

pix^") = pix") . (6.11) 
Expanding both sides of Eq. (6.11) using Eq. (2.1), we get 

p{x'')^Po{x^ + eSp^ix^) + 0{e^) , (6.12) 

= Po(x°) + e {-p'o{x')eAxn + S^iixn) + 0{e') . (6.13) 

Note that we use the same background solution po(^) in both expressions. 
Thus we obtain the transformation rule at first order 

Spi = 5p, + p',ei. (6.14) 
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Another invariant is the hne element ds"^, which allows us to deduce the trans- 
formation properties of the metric tensor by exploiting the invariance of ds^, 
i.e., 

ds'^ = g^ydx^dx^ = g^^dx^dx" , (6.15) 
We here will not follow this approach further here, but see e.g. [66,104]. 

6.2 Four-scalars 

We now return to the active approach by studying the simplest tensorial quan- 
tity, the four-scalar. Examples of four-scalar are the energy density, p, and the 
scalar field ip, and we shall use the former below. 

Prom Eqs. (2.1) and (2.2) we immediately get the perturbed four-scalar up to 
second order ^ 

Po + 5pi + -5p2, (6.16) 
using the energy density as an example. 

6.2.1 First order 

Before we can study the transformation behaviour of the perturbations at first 
order, we split the generating vector into a scalar temporal part cti and a 
spatial scalar and vector part, /3i and 7^*, according to 

er=(ai,/3i/ + 7iO , (6.17) 
where the vector part is divergence-free dk^^^ — 0. 

Under a first-order transformation a four scalar, here the energy density, p, 
then transforms from Eqs. (6.5) and (A. 8) as, 

5pi^ 5pi + p'^ai. (6.18) 

The first-order density perturbation is fully specified by prescribing the first 
order temporal gauge or time slicing, a^. 

6.2.2 Second order 

At second order, as at first order, we split the generating vector ,^2 i^^^o ^ 
scalar time and scalar and vector spatial part, similarly as at first order, as 

e2 - p^: + 1^) , (6.19) 



21 



where the vector part is divergence-free dk')-^ — 0. We then find from Eqs. (6.6) 
and (A. 8) that a four scalar transforms as 



+ (2<5pi+pW,,(/3i,'= + 7i1- 



(6.20) 



We see here aheady the couphng between vector and scalar perturbations in 
the last term through the gradient and '-y^. The gauge is only specified once 
the scalar temporal gauge perturbations at first and second order, ai and 0:2, 
and the first order spatial gauge perturbations, /3i and 7^', are specified. 

6.3 Four-vectors 

We now turn to four- vectors and their transformation properties. Of particular 
interest in cosmology is the unit four-velocity which we defined in Section 
4 above ^ . 

6.3.1 First order 

A four- vector transforms at first order, using Eqs. (6.5) and (A. 9), as 



where we used the fact that in a FRW spacetime background quantities are 
time dependent only. 

For the specific example of the four- velocity, defined in Eq. (4.3), we find. 



Using the transformation of the metric perturbation Si,, given below in Eq. (6.30), 
and using the decompositions given in Section 2.2.1 above, we get the trans- 
formations for the scalar and vector parts, respectively, at first order 



Note that under a gauge transformation the vector field normal to the constant- 
T] hypersurface, defined in section 3.1, will be replaced by a new vector field 
normal to the constant-?? hypersurfaces, therefore it is not particularly helpful to 
study the transformation of n^. We will study the gauge transformation of the 
metric tensor in section 6.4, from which we can derive the transformation rules for 
the hypersurface-orthogonal field and its expansion, shear, etc. 




vii + Bii = vii + Bii - ai^i . 



(6.22) 
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^i = v^-l3[, (6.23) 
5^ = <ec(i)-7/, (6.24) 

and the perturbed temporal part of does indeed transform as a scalar. 



6.3.2 Second order 

At second order we find that a four-vector transforms, using Eqs. (6.6) and 
(A.9), as 



5^2/. = 5^2^ + W('o)^a2 + W(o)o«2,M + W(o)m"? + (6-25) 
+2^'o)oaiai,^ + W(o)o (Ci «i,mA + "i.aCi, J + 2 {^U^^Ai + SU^x^t) ' 

where as before we used for the background U(o)^ = U(o)ti{v) ^-^d U(o)i = 0. 

Focusing again on the four- velocity, Eq. (4.3), and following a similar pro- 
cedure as at first order, we find that the second order combined scalar and 
vector spatial part transforms as 



V2i = V2i - + ^yi > (6-26) 

where contains the terms quadratic in the first order perturbations and is 
given by 



= (201 + a[ + 2Ha,) - a,Cu 

-^i^u,k + ^lUk - 2ai {v[, + Hvu) + 2vu,kei - , (6-27) 

and we already substituted for the transformation of the metric perturbation 
B2i, given below in Eq. (6.48). Decomposing then the second order velocity 
transformation, Eq. (6.26), into scalar and vector parts, we get the transfor- 
mations as 



"^vec(2) 



(6.28) 

(6.29) 
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6-4 Tensors 



6.4-1 First order coordinate transformation 

We can now calculate how the first order metric perturbations change under 
a gauge transformation. We get the change of the 6^^^ goo and hence the lapse 
function 0i immediately from Eqs. (6.5) and (A. 10), since this component of 
the metric is scalar in nature. The change of the S^^^goi slightly more involved, 
since this component contains scalar and vector perturbations. We therefore 
have to compute the overall transformation of this metric component using 
Eqs. (6.5) and (A. 10), and then split the result for into its scalar part, Bi, 
and its divergence-free part —Si. We therefore get for the combined part Bn, 

B,, = B,i + ^[,-ai,i, (6.30) 
and taking the divergence gives for the scalar part, 

V^Bi^V^Bi + V^(3[-V^au (6.31) 

which, after "removal" of the Laplacian gives the transformation behaviour of 
Bi. We can then subtract the scalar part from Eq. (6.31) and are left with the 
vector part. The results are given below. 

To get the change of the metric functions in the spatial part of the metric 
under a gauge transformation, we again first use Eqs. (6.5) and (A. 10) to get 
transformation of the spatial part of the metric S^^^Qij, or Cuj, 

2Cuj = 2Cuj + 2HaiSij + + ^ij^i , (6.32) 
where we reproduce Eq. (2.12) above for convenience, 

2Ciij — —2i>iSij + 2Ei^ij + 2Fi(jj) + huj . 
Taking the trace of Eq. (6.32) we get 

-3VJi + V^Ei = -3V^i + V^^^i + SHai + V^A . (6.33) 

Now applying the operator d^d^ to Eq. (6.32) we get a second equation relating 
the scalar perturbation ipi and Ei, 

-3V2Vi + V'V'^i = -3VVi + V'^V^Ei + mVai + V . (6.34) 
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Taking the divergence of Eq. (6.32) we get 



2C^^^ = 2Ci,^./ + 2na,,i + V^^H + • (6.35) 

Substituting in our results for and Ei we arrive at 

V'Fh = V'Fh + VV- (6.36) 



We can sum up the transformations of the first order metric perturbations we 
have from the above, first for the scalars as 



01 = 01 + Ti-Oii + a\ , (6.37) 

i^i — i^i — Ti-O-i , (6.38) 

Bi^Bi-ai + P[, (6.39) 

Ei^Ei+Pi, (6.40) 

and for the vector perturbations as 

= (6.41) 
Fi^ = F/ + 7/. (6.42) 

The first order tensor perturbation is found to be gauge- invariant, 

hiij = hiij . (6.43) 
by substituting Eqs. (6.37) to (6.42) into Eq. (6.32). 

For later use, we note that the scalar shear potential, ai = E[ — Bi, defined in 
Eq. (3.15) above, and the combination Vi+Bi corresponding to the momentum 
scalar, transform as 

51 = (7i + ai , (6.44) 
n + Bi^vi + Bi-ai. (6.45) 

6.4-2 Second order gauge transformations 

The metric tensor transforms at second order, from Eqs. (6.6) and (A. 10), as 
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^„(2) _ r„(2) , (0) c\ , , „(0) t 



.(0) CA 



(0)cA 



-,(0)cA 



(l)cA 



, „(0) tAta , (0) tA ta , oLW <r«<rA , (0) ta^A , ^WtA 

{ii + .) + gf} {e, + ei jt ,) • (6.46) 



As at first order, in the previous subsection, we get the transformation be- 
haviour for the second order lapse function 02 straight from the — compo- 
nent of Eq. (6.46), which gives 



02 = 02 + 'Ha2 + a2 + 0L\ 



+2ai (a/ + 2di) + (a/ + Hai + 20i 



(6.47) 



The combined scalar and vector — i metric part transforms from Eqs. (6.6) 
and (A. 10) as 



B2i = B2i + ^2i ~ (^2,i + ^Bi , 



(6.48) 



where i?2i is similarly vector and scalar combined, and we defined to 
contain the terms quadratic in the first order perturbations, as 



To get the transformation behaviour of the vector and the scalar part sepa- 
rately, we take the divergence of Eq. (6.48) and find after applying the inverse 
Laplacian, the transformation scalar part S2, 



E2 = 52-a2 + /?2 + V-'A'B\, 



(6.50) 



or explicitly 
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Ic k T-72 : 



. (Ci^ k - < k) - <^i,k^%' - <ki' , j ■ 

The vector part is then simply found by subtracting the scalar part from 
Eq. (6.48), and is given by 



(6.51) 



821 — 821 — 72* — + V ^X-B^M ■ 



(6.52) 



We now turn to the transformation behaviour of the perturbations in the 
spatial part of the metric tensor. We can follow here along similar lines as in 
the linear case. However, the task is made more complicated not only by the 
size of the expressions but more importantly by the fact that now we will have 
to let inverse gradients operate on products of first order quantities. 

The perturbed spatial part of the metric, C2ij, transforms at second order as 



2(^211 — 2C2ij + 2Ha2Sij + ^21,3 + ^21,1 + ^ij , 



(6.53) 



where we defined Xij to contain the terms quadratic in the first order pertur- 
bations as 



A'jj = 2 



n'^ + —jaj + n + a^^k^f 

Hu,kCi' J + iij,kii' , + iuaij + ■ (6.54) 



+ 2 {Buaij + Bijai^i) 



The perturbed spatial part of the metric, C2ij, is decomposed in Eq. (2.12) 
above into scalar, vector, and tensor part, which we reproduce here at second 
order, 

2C2ij = —2ip2Sij + 2£^2,ii + 2F2(ij) + /i2ij ■ 
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Taking the trace of Eq. (6.53) we get 



-3ilj2 + V^E2 = -2>iJ2 + V^E2 + ma2 + ^132 + -X\ , 



where we find X\ to be 



(6.55) 



+2 



(6.56) 



Now applying the operator d'^d^ to Eq. (6.53) we get a second equation relating 
the scalar perturbations ■02 and £'2, 

- V V2 + W^W^E2 = -W^iJ2 + W^W^E2 + 7t:v2a2 + W2/52 + ^ A'^^' , (6.57) 

This gives the second-order scalar metric perturbations 

1. 



and 



= v^2 - na2 - -x\ + , 



E2^E2 + P2 + ^V-'VX'\,^ - \V-^X\ 



Taking the divergence of Eq. (6.53) we get 

2C2,,.^' = 2C2,,-^' + 2na2, + V^e2^ + V2/52,, + X,^^ 
Substituting in our results for ■02 and E2 we then arrive at 

V'F2. = V^F2, + + A",,, ^ - V-^A-^',,,, . 
Finally we obtain the second-order vector metric perturbation 

hi = F2i + 72i + V-^A",,, - V-^ V-^A''=',,,, . 



(6.58) 
(6.59) 
(6.60) 

(6.61) 
(6.62) 



We can now turn to the tensor perturbation at second order. Substituting our 
previous results for ■02, -£'2, and ¥21 into Eq. (6.53) we get 



h2ij — h2ij 



-^--^X^ij - [X^k, ^ j + ^jk, ^ i) ■ 



,klij 



(6.63) 



28 



Although the second-order tensor transformation h2ij is not dependent on the 
second-order part of the gauge transformation, does depend on first 

order quantities quadratically. The tensor metric perturbations are no longer 
gauge-invariant at second and higher order. 

6.4-3 The large scale or small k limit 

Prom Eqs. (6.20), (6.58), and (6.47) we see that on large scales, where gradient 
terms can be neglected, the definition of the second order perturbations in 
the "new" coordinate is independent of the spatial coordinate choice (the 
"threading") at second order in the gradients. It is therefore sufficient on 
large scales (at 0{k'^ j) to specify the time slicing by prescribing cti and 0:2, in 
order to define gauge- invariant variables [110,98,107]. The procedure to neglect 
the gradient terms, is explained in detail in Ref. [98]. For the approximation 
to hold one assumes that each quantity can be treated as smooth on some 
sufficiently large scale. Formally one multiplies each spatial gradient di by a 
fictitious parameter k, and expands the exact equations as a power series in 
k, keeping only the zero- and first-order terms, finally setting k — 1. 



7 Gauge- invariant veiriables 

The notion of invariance under coordinate reparametrisation is central to Ein- 
stein's theory of general relativity. This is both a blessing and a curse in the 
study of cosmological perturbations. We are free to pick coordinate systems 
best adapted to the problem at hand, but we also obtain apparently different 
results depending upon this arbitrary choice of coordinates. This is the gauge 
problem. 

Ultimately physical observables are not dependent on the choice of coordinate, 
though observables may be different for different observers. All one can do is 
to specify quantities unambiguously, such that they have a gauge-invariant 
definition. This is not the same as gauge independence. A quantity like the 
tensor metric perturbation, hnj, is truly gauge independent at first order in 
that the tensor part of the metric perturbation is the same in all gauges. The 
scalar curvature perturbation, ipi, on the other hand is intrinsically gauge- 
dependent. It is different under different time slicings. (Indeed in the spatially 
fiat slicing the curvature perturbation is zero by construction.) One can con- 
struct gauge-invariant combinations, which may be referred to as the gauge- 
invariant curvature perturbation, but they only correspond with the curvature 
perturbation in one particular gauge. As a result one can find in the litera- 
ture many different gauge-invariant curvature perturbations corresponding to 
the many different choices of gauge, such as ^, C and TZ, corresponding to 
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the curvature perturbation in the longitudinal, the uniform density, and the 
comoving gauge, respectively, to name just three. 

In this section we shall show how different gauge invariant combinations of 
otherwise gauge dependent quantities can be constructed by fixing the other- 
wise arbitrary coordinate transformations at first order and beyond, yielding 
gauge invariant definitions of the physical perturbations in specified gauges. 
Residual gauge degrees of freedom only remain in cases where the coordinate 
choice is not unambiguously fixed, as in the synchronous gauge. 

At first order the tensor metric perturbation, huj, is not affected by the map- 
ping, or by the change of coordinate system, and hence is gauge-invariant. 
Thus we only have to construct gauge- invariant scalar and vector perturba- 
tions at first order. However at second order the tensor part of the metric 
perturbation also becomes gauge-dependent. 

At first order we can define scalar and vector type gauge-invariant variables 
independently of each other, but matters are more complicated at second 
order. Whereas we can still specify the "proper" second order scalar and vector 
slicings and threadings independently, we now also have to specify the first 
order gauge functions of both types simultaneously. 

We could specify different gauges at first and second order, but we would 
be losing the physical interpretation of the quantity thus constructed. We 
therefore choose the same gauge at first and second order, and at first order 
the same physical gauge condition for the scalar and vector gauge functions. 
It is however sometimes necessary to combine different temporal and spatial 
gauge conditions. For example imposing the uniform density condition only 
specifies the slicing, and we are free to combine it with a fiat threading. 



7. 1 Longitudinal gauge 
7.1.1 First order 

The gauge-dependence of the metric perturbations lead Bardcen [7] to pro- 
pose that only quantities that are explicitly invariant under gauge transforma- 
tions should be considered. By studying the transformation Eqs. (6.37-6.40), 
Bardeen ^ constructed two such quantities [7] 



^ In Bardeen's notation these gauge-invariant perturbations are given as ^> = 
^aQ^^^ and * = -^hQ^^^- 
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$ = 01 + n{B, - E[) + (Si - E[y , (7.1) 

^ = 'ijji-n(B,-E[) . (7.2) 

These turn out to coincide with the scalar metric perturbations in a partic- 
ular gauge, called variously the orthogonal zero-shear [7,66] , conformal New- 
tonian [21,101] or longitudinal gauge [123]. It may therefore appear that this 
gauge is somehow preferred over other gauge choices. However any unam- 
biguous choice of time-slicing and threading can be used to define exphcitly 
gauge-invariant perturbations. The longitudinal gauge of Ref. [123] provides 
but one example. 

The two scalar gauge functions, a and f3 defined in Eq. (6.17), which represent 
different choices of time-slicings (choice of spatial hypersurfaces) and threading 
(choice of spatial coordinates on these hypersurfaces) respectively, allow two 
of the scalar metric perturbations to be eliminated, implying that the two 
remaining gauge-invariant combinations should then be gauge-invariant. If we 
choose to work on spatial hypersurfaces with vanishing shear, we find from 
Eqs. (6. 39), (6. 40) and (3.15) that the shear scalar transforms as oi = cxi -|- cti 
and this implies that to obtain perturbations in the longitudinal gauge starting 
from arbitrary coordinates we should perform a transformation 

au = -(Ti =Bi-E[. (7.3) 

This is sufficient to determine the geometrical perturbations 0i, V"!, en or other 
scalar quantities on these hypersurfaces. In addition, the longitudinal gauge 
is completely determined by the spatial gauge choice Eu = [and hence from 
Eq. (7.3) Bu = 0] which requires from Eq. (6.40) 

Pu = -El . (7.4) 

The remaining scalar metric perturbations, (pi and ipi, are given from Eqs. (6.37) 
and (6.38) as 



<P^ = <Pi + n{Bi - E[) + {Bi - E[)' , (7.5) 
il^u = i^i-n{Bi-E[) , (7.6) 

Note that and are then identical to $ and ^ defined in Eqs. (7.1) 
and (7.2). 

The fluid density perturbation, and scalar velocity, Vi, are given from 
(6.18) and (6.23) 



5pu^5pi + p'o{Bi-E[) , (7.7) 
€{i^vi + E[. (7.8) 
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These gauge-invariant quantities are simply a gauge-invariant definition of the 
perturbations in the longitudinal gauge. 

This gauge is widely used, for example, throughout Ref. [123]. It has also 
proven useful for calculations on small scales, since it gives evolution equations 
closest to the Newtonian ones, e.g. Ref. [54]. Recently it has also become 
popular in backreaction studies, e.g. [19,71,170]. After imposing the gauge- 
conditions the metric tensor is diagonal, which simplifies many calculations, for 
example the derivation of the governing equations of the Boltzmann-hierarchy. 
Moreover we shall show in Section 8 that in many cases of physical interest 
(in the absence of anisotropic stress) one finds $ = and there is only one 
variable required to describe all scalar metric perturbations. 

However, it can be difficult to define quantities in this gauge in the super- 
horizon limit, since in the super- horizon limit the shear vanishes and hence 
numerical instabilities can occur on large scales in the longitudinal gauge, see 
e.g. Refs. [101,60,82]. 

The extension to include vector and tensor metric perturbations is called the 
Poisson gauge [21,26]. Tensor metric perturbations are automatically gauge 
independent at first order (and hence gauge invariant). Eliminating the spa- 
tial part of the contravariant vector field in Eq. (3.5) requires both Bi 

and Se = which from Eq. (6.41) fixes the vector part of the spatial gauge 
transformation (6.17) 

Yu = J Sldv + Cl{x=), (7.9) 

up to an arbitrary constant 3- vector €{ which depends on the choice of spatial 
coordinates on an initial hypersurface. The remaining vector metric perturba- 
tion is 

Fu = Fi' + I Sidrj + Ci{x^) . (7.10) 

7.1.2 Second order 

It is possible to extend the longitudinal, or Poisson, gauge to higher orders. 
The principle for constructing gauge invariant variables remains the same. We 
use a physical choice of gauge to specify the vector field generating the 
transformation (6.7) from an arbitrary gauge [110,124]. Requiring first that 
— fixes the scalar part of the spatial gauge using Eq. (6.59), which gives 

P2e = -E2 - + \v-'X\ . (7.11) 

Note that having already imposed the Poisson gauge at first order, ai, j3i 
and 7ij are fixed by Eqs. (7.3), (7.4) and (7.9), and thus so is A'j^, given in 
Eq. (6.54). 
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Requiring that the scalar part of the perturbed shift function vanishes, — 
0, then sets the temporal gauge, 0:2^ using Eq. (6.50), while requiring that the 
vector part vanishes, F2 = 0, can be used along with Eq. (6.62) to fix the 
vector part of the spatial gauge, 72^, up to a constant of integration as at first 
order, Eq. (7.9). 

We then obtain gauge invariant definition of and other perturbations at 
second order by substituting these specific gauge transformations into Eq. (6.47) 
and (6.58) to obtain 



+2au (au + 20i) + ^uk {au + Hau + 20i) 



Uk 



a 



u. 



— 2Bik — ^ 



u 



i^2e = V'2 - ^-012 



1, 

4 



(7.12) 
(7.13) 



where Xuj is denotes the quadratic first order terms in Eq. (6.54) using the 
longitudinal gauge transforms au and ^^j. 



The tensor (transverse, traccfrcc) part of the metric perturbation at second 
order, hij in Eq. (6.63), is not afi^ected by the second order gauge transforma- 
tions a2 and ^2i, but it does depend on the choice of gauge at first order, ai 
and ^ij. Thus we need to include the corresponding first-order gauge defini- 
tion to obtain a gauge invariant definition of the tensor metric perturbations 
at second order. 



In particular, recent work on the generation of gravitational waves at second 
order [117,126,5,18,16] has calculated the resulting tensor mode in the Poisson 
gauge. To give a gauge invariant definition of the tensor perturbation in the 
Poisson gauge one needs to explicitly include the transverse and tracefree (ten- 
sor) part of the second order gauge transformation from an arbitrary gauge. 
From Eqs. (6.63) we obtain the gauge invariant definition of the tensor metric 
perturbation in the Poisson gauge: 



1, 
2 



+-V '^X^k,ij - ^ ^ i^iik, 



ejk, 



0- 



(7.14) 
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1.2 Spatially flat gauge 



7.2.1 First order 

An alternative gauge choice, defined purely by local metric quantities is the 
spatially flat or uniform curvature gauge [66,62,63,65,161], also called the off- 
diagonal gauge [27]. In this gauge one selects spatial hypersurfaces on which 
the induced 3-metric on spatial hypersurfaces is left unperturbed by scalar 
or vector perturbations, which requires ipi — Ei — and F^ = 0. Using 
Eqs. (6.38), (6.40) and (6.42) this corresponds to a gauge transformation (6.17) 
where 

«lflat = ^, Aflat = -^1, 7iflat = -^l'- (7-15) 

The gauge-invariant deflnitions of the remaining scalar metric degrees of free- 
dom are then from Eqs. (6.37) and (6.39) 



01flat = 01+^l+(^) , (7.16) 

B^,^B,-E[-'^. (7.17) 

At 

These gauge-invariant combinations were denoted A and B by Kodama and 
Sasaki [66]. The gauge- invariant deflnition of the remaining vector metric per- 
turbation is the time derivative of the vector metric perturbation in the Poisson 
gauge (7.10): 

S^i^Sl + Fl' = K/. (7.18) 



Perturbations of scalar quantities in this gauge, such as the density perturba- 
tion, have gauge invariant definitions from Eq. (6.18): 

Spm^t^ Spi + p'o^ , (7.19) 

and the velocity potential (6.23) is given from Eq. (6.23): 

vl^^t ^Vi + E[. (7.20) 



The shear perturbation in the spatially fiat gauge is given by Oia^ = — -Bisat- 
This is closely related to the curvature perturbation in the zero-shear (longi- 
tudinal) gauge, ipi — ^, given in Eqs. (7.2) or (7.6), 

i>u 

-Biflat = = ■ (7-21) 
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Gauge-invariant quantities, such as Smat or ipi£ are proportional to the dis- 
placement between two different choices of spatial hypersurface, 



-Blflat = — 



n 



QJlflat — Olu , 



(7.22) 



which would vanish for a homogeneous cosmology. 



In some circumstances it is more convenient to use the spatially-flat gauge- 
invariant variables instead of those in the longitudinal gauge. For instance, 
when calculating the evolution of perturbations during a collapsing "pre Big 
Bang" era the perturbations 0iflat and -Biflat may remain small even when $ 
and ^ become large [27,39] . On the other hand the metric perturbation -BiAat 
grows on large scales in radiation or matter dominated eras, Siflat oc Ti"^ oc 77, 
when the longitudinal gauge metric perturbation "if remains constant. 



Note that the scalar field perturbation on spatially fiat hypersurfaces. 



(7.23) 



is the gauge-invariant Sasaki-Mukhanov variable [145,120], often denoted Q. 
7.2.2 Second order 

At second order we get from the gauge condition ip2 — using Eq. (6.58) 



1 



Q;2flat 



Y7-2 yij yk 
^ ^flat,ii ^flatfe 



(7.24) 



where we get Aflatij from Eq. (6.54) using the first order gauge generators given 
above, as 



'^flatii — 2 



lij 



+4CHj,fe^fflat + (^Ciife + 6flati,fe) ClflatJ + i^Cijk + 6flatj,fe) Clflat,? 



^Iflati) 



where we define 



+ (Clflat(iJ) + 4^Clflat(i,j)) + 2CmatClflat(i,j)fe + 2Clflatfc,iC 



j Iflati 



k 

Iflat ,j ' 

(7.25) 
(7.26) 
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The trace of Eq. (7.25) is then 



n 



+4Cf + 4 (2C^' + ek,, ') - 2V^£;i,fc^ffl,, (7.27) 



n 



7i 



As an example of a second-order scalar metric perturbation we give the lapse 
function in the flat gauge 



^2flat 



n 



^2+ (^-^) ^2 



(7.28) 



Y7-2 ykl/ ykl j_ ( %/ ^ 

^ ^flat,fei ^flatfe "Tin. — 



4 ^ 



+ 



+ 



,k 



7i 
- 2'HBik 



■Slflat 



»lflat 



The second order tensor perturbation is in the flat gauge 
~ 1 1 

+ -V"^A'4tfc,ij - ('^flatffcj + '^flatjTci) ■ (7.29) 

As an example of a matter variable we choose the energy density, which in the 
flat gauge is 



<^P2flat = + ^^2 + ^ - ^Lk) 



+ 



PoVi + Po(V'l-^V'il+2Wi 



(7.30) 

+ (2<5pi + ^Vi) 

,k 



To show the relation between gauge-invariant perturbations defined in the flat 
gauge and those previously defined in the longitudinal gauge we note that in 
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the longitudinal gauge we have 



2C 



Uij 



2^% + 2F, 



U{i,j) -t- hiij , 



and the first-order gauge shifts to the flat gauge, Eqs. (7.15), become 



QilflatU 



n 



IflatI 



0, 



at 



(7.31) 



(7.32) 



Substituting these into Eq.(7.25) we obtain 



n 



\ k k 

2^^kOij — 2Fi£(jj);j — hij^k ) Fii — 4hk{iFuj) 

2 / k k 

~y^^Fu(ij) + 2Fu Fu(ij)k + '^Fuk,iFuj ■ 



(7.33) 



If we can neglect the first-order vector and tensor perturbations, Fm and huj, 
then we have 



(7.34) 



Similarly, neglecting first vector and tensor perturbations, i.e. setting F^ — 
Sii = hiij = and hence only considering scalar perturbations, the energy 
density at second order, Eq. (7.30), simplifies to 



dp2flat = 0P2 + —^2 + ^ 

-2 [Sp, + ^ij?j ^E,>: 

+Ai^E,,kiE,':' + V'E^,kE,'; - 2^iV^ + 2Er 



{pI + 2np'^) ^1 + p'^ (2i^[ - ^^i] + 2 Wi 



I ij 



n 



n 



2^1 1 ^ + 2^ir' ) - E^f'E,,^ - E,]^E,l + i^[2B, + E[- 



n 



n 



This expression was first derived in Ref. [110], however with a different, in- 
correct sub-horizon part [133]. On super-horizon scales, where gradient terms 
can be neglected, we recover the expressions given in Ref. [105]. 
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1.3 Synchronous gauge 



The synchronous gauge is defined by ^ — Bi — 0, so that the proper time for 
observers at fixed spatial coordinates coincides with cosmic time in the FRW 
background, i.e., dr = adrj. This simphfies dynamical equations as the time 
derivatives can be directly related to proper time derivatives. This gauge is 
very popular for numerical studies, and used in many Boltzmann solvers such 
as CMBFAST [154]. It is also popular in the older literature [177,131]. 

The gauge condition at first order is 0i = Bi^ — 0, which from Eqs. (6.37) 
and (6.39) gives 



aisyn = - ^ (/ «<^i^^ - Ci{x')^ , (7.36) 
/3isyn = / (aisyn " ^i) dr} + Ci{x') , (7.37) 
lly^^Jsidrj + Ciix'). (7.38) 

This docs not determine the time-slicing unambiguously and we are left with 
two arbitrary scalar functions of the spatial coordinates, Ci and Ci. Note that 
+ Cii afi:ects only the labelling of the coordinates on the initial spatial 
hypersurface, but Ci affects scalar perturbations on spatial hyp crsurf aces. We 
are left with two non-zero geometrical scalar perturbations, 



V'lsyn = ^1 + 7 (/ ''^idr] - C{x')j , (7.39) 

O^lsyn = CTi + ftigyn - Bi , (7.40) 



and the matter variables are 



<^Pisyn ^5pi-^(^J a<Pidr) - C{x')^ , (7.41) 

^^l^^n ^Vi + Bi- aisyn ■ (7.42) 

Thus it is not possible to define gauge-invariant quantities in general using 
this gauge condition [114]. This gauge was originally used by Lifshitz in his 
pioneering work on perturbations in a FRW spacetime [87] (see also Ref. [75]). 
He dealt with the residual gauge freedom by eliminating the unphysical gauge 
modes through symmetry arguments. 

To remove the ambiguity, we can follow Ref. [28] and choose the initial velocity 
of cold dark matter to be zero, Vi^dm. = 0, which fixes the residual gauge 
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freedom 

Ci{x) = a(?;icdm + Bi) . (7.43) 
Note that for pressureless matter, momentum conservation equation ensures 
that a(f icdm + Bi) is a constant (see Section 8). 



7.4 Comoving orthogonal gauge 



The comoving gauge is defined by choosing spatial coordinates such that the 
3- velocity of the fluid vanishes, Vi = 0. Orthogonality of the constant-?] hyper- 
surfaces to the 4- velocity, u'^, then requires Vi + = 0, which shows that the 
momentum vanishes as well. Prom Eqs. (6.39) and (6.23) this implies 

ttlcom = Vi + Bi , 

Acorn = / V,d7J+C{x'), (7.44) 

where C{x^) represents a residual gauge freedom, corresponding to a constant 
shift of the spatial coordinates. All the 3-scalars like curvature, expansion, 
acceleration and shear are independent of C(a;*). Applying the above transfor- 
mation from arbitrary coordinates, the scalar perturbations in the comoving 
orthogonal gauge can be written as 



Alcorn = <Pi + n{vi + B,) + {v[ + B[) , (7.45) 
n = i^^ = iJi-n{vi + Bi) , (7.46) 
S^^vi + E[. (7.47) 

Defined in this way, these combinations are gauge-invariant under transfor- 
mations of their component parts in exactly the same way as, for instance, $ 
and ^' defined in Eqs. (7.1) and (7.2). 

Note that the curvature perturbation in the comoving gauge given above, 
Eq. (7.46) was used (with a constant pre-factor) by Lukash [91]. It was later 
employed by Lyth and denoted TZ in his seminal paper, [92], and in many 
subsequent works, e.g. [83] and [86]. 

The density perturbation on the comoving orthogonal hypersurfaces is given 
by Eqs. (6.18) and (7.44) in gauge-invariant form as 

Spicom = Spi + po (vi + 5i) , (7.48) 

and corresponds to the gauge-invariant density perturbation emEoQ'^'^^ in the 
notation of Bardeen [7]. The gauge- invariant scalar density perturbation A 
introduced in Refs. [25,46] corresponds to Spi^j^^^/po. 
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If wc wish to write these gauge-invariant qiiantities in terms of the metric 
perturbations rather than the velocity potential then we can use the Einstein 
equations, presented in section 8, to obtain 

(7.49) 

In particular we note that we can write the comoving curvature perturbation, 
given in Eq. (7.46), in terms of the longitudinal gauge- invariant quantities as 

which coincides the quantity denoted ^ by Mukhanov, Feldman and Branden- 
berger in Ref . [123] . 

By comparing the definitions of the energy-momentum tensor for a single fiuid 
and a single scalar field in Sections 4.1 and 5.1 we can relate the Vi + Bi to 
6(pi, which allows one to rewrite the definition of the comoving curvature 
perturbations, Eq. (7.46), as 

7^ = V'l + ^Sipi . (7.51) 

From the definition above we immediately that the comoving curvature per- 
turbation is related to the field fluctuation on flat slices, deflned in Eq. (7.23), 
by 

S^t = ^T^- (7-52) 
For extensions to the multi-field case see Section 8.2.4. 



7. 5 Total matter gauge 

This gauge is also known as the velocity orthogonal isotropic gauge [66] but 
here we follow the terminology of Ref. [84] . It is closely related to comoving 
orthogonal and longitudinal gauges. 

To fix the temporal gauge we require the total momentum potential on spatial 
hypersurfaces to vanish 

v{ + Bi^O. (7.53) 

In addition we require Ei — and Fij = 0, which fixes the spatial gauge. 
These require 

aitom ^Vi + Bi, /3itom = -El , 7itom = "^i' ■ (7-54) 
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We therefore get the metric perturbations in the total matter gauge related 
to the comoving orthogonal and longitudinal gauge perturbations: 



^Itom — 01com ) 
V'ltom — V'lcom , 
-Sltom = —Vi£ , 



(7.55) 
(7.56) 
(7.57) 



and for the matter quantities we get in the total matter gauge 



^Pltom — (^Plcom ) 



t'ltom = Vu ■ 



(7.58) 

(7.59) 



Note that in the total matter gauge velocity potential is not identically zero 
(unlike in the comoving orthogonal gauge), but equal to the shear potential, 

vuma = ^uma = — -Bitom, which also coiucides with velocity potential in longi- 
tudinal gauge. 

7.6 Uniform density gauge 

Alternatively we can use the matter to pick out a foliation of uniform density 
hypersurfaces on which to define perturbed quantities. 

7.6.1 First order 

Using Eq. (6.18) we see that 5pi — Q implies a temporal gauge transformation 



On these hypersurfaces the gauge-invariant curvature perturbation is [40,114] 



The sign is chosen to coincide with C, defined in Refs. [9,8] ^ . There is still the 
freedom to choose the spatial gauge. In particular we can choose either B, E 
or V to be zero and thus fix (5. 



^ Note, that Ci defined in Eq. (7.61) is related to the curvature perturbation Csbb 
defined in Ref. [143] by Csbb = 3Ci. 



5pi 



(7.60) 





(7.61) 
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Note that (i is simply related to the density perturbation in the flat gauge, 
Eq. (7.19) by 

-Ci = ^ Wt ■ (7.62) 

Po 

The curvature perturbation in the uniform-density gauge is also closely related 
to the comoving curvature perturbation (7.46). At first order we have 

Ci = -T^i r^Picom , (7.63) 

Po 

where 5picom is the comoving density perturbation (7.48). In Section 8 we shall 
use the Einstein equations to show that and TZi differ only by an overall 
minus sign in the large scale limit where the comoving density perturbation 
vanishes. 



7.6.2 Second order 

The transformation behaviour of scalars at second order, Eq. (6.20), allows us 
to define the temporal gauge corresponding to uniform density hypersurfaces 
as 

a25p = - ^ (p'o'ai + p'^a[ + 26 p[) - \ {26 p, + p'^a,) ^ (p,^ ^ + 7i 

fin fli-i fin ^ 



Po Po 



Po 



(7.64) 

Using then the definition of uniform density hypersurfaces at first order, 
Eq. (7.60), and choosing a spatially flat threading by using Eq. (7.15), we 
flnally get 



a25p - -7 
Po 



Po ^ ' . 



(7.65) 



Using Eq. (6.54) we flnd 



ijSp 



n (i + 3ci) j _ ^Asp, + 6p,,,i,' 



(7.66) 



+4 



6p 



~r iplij + 27iCijj) + Ciij^kil' + (4Ciifc + Clflati,fe) Clflat J 



Po 



+ (4Cijfe + Clflatj.fe) Clflat.i - -J ^Pl,i (2-Blj + Clflatj) + ^Phj + Clflati) 



Po 



S p ^ 

-2— (Clflat(jJ) + 47YClflat(i,j)) ^Spi,i6pij + 2^iflat6flat(jj)fc + 2Clflatfe,iC 

Po Po 



k 

IflatJ ' 
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where we choose a flat threading by defining, 



(7.67) 



Prom Eq. (7.66) we get for the trace in the uniform density gauge 



n 

7o 



\ Po / Po 



+4 



Po 



+4 (2Cf + em.t,') Ciflat(.,o - - (2^1' + ^m.tk) Spi,k 

Po 



Po 



Po 



(7.68) 



Finally, from Eq. (6.58) we get for ^2, the curvature perturbation on uniform 
density hyp ersurf aces. 



-C2 = -02 + — 

Po 



8P2 - ^^Spi + 

Po 



k^,. + 7V-''^^.-(7-69) 



The second order tensor perturbation is in the uniform density gauge 



h25pij — h2ij + Xspij + 2 ^'^Sp,kl - '^Spk) ^ij + 2^ ^'^Sp,klij 

+-V"^A'/pfc_jj - V"^ {^spik, j + '^Spjk, i) ■ (7.70) 



The curvature perturbation on uniform density hypersurfaccs, defined in Eq. (7.69), 
simplifies considerably if we neglect first order vector and tensor perturbations, 
i.e. setting = Su = huj = and hence only considering scalar perturba- 
tions. 
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-C2 = V2 + —Op2 Y 

Po Po 



Po 



Po 



(7.71) 



Pb 



Pb 



Pb 



\ Po / Po \ Po / 1, 

This expression was first derived in Ref. [110], however with a difi^erent, in- 
correct sub-horizon part [133]. On super-horizon scales, where gradient terms 
can be neglected, we recover the expressions given in Ref. [105]. 



8 Dynamics 



In this section wc give the Einstein equations governing the evolution of the 
FRW background and perturbations in general relativity. This will allow us to 
derive some key properties of the perturbation variables, such as the conserva- 
tion of the curvature perturbation C, on super-horizon scales in the adiabatic 
case. 

In general relativity the Einstein equations relate the local spacetime curvature 
to the local energy-momentum: 

G^, ^SttCT;,. (8.1) 

In more general theories of gravity we can still equate the local spacetime 
curvature, G^,/, with an effective energy- momentum, though this may not be 
simply related to the energy-momentum tensor derived, say, from the matter 
Lagrangian. Moreover, many modified gravity theories, including Brans-Dicke 
gravity or higher-order theories, may be rewritten in terms of general relativity 
plus non-minimally coupled matter fields through a conformal rescaling of 
coordinates [102,173]. In this review we will restrict our analysis to general 
relativity. 

We can project the tensor equation (8.1) into components tangent to and 
orthogonal to the timelike 4- vector field, defined in Eq. (3.4), which defines 
the coordinate system (see section 3). This gives two constraint equations for 
the metric perturbations, which we will refer to as the energy and momentum 
constraint equations. We also have two evolution equations driven by the trace 
and trace-free parts of the pressure. Through the Bianchi identities, V^G^ = 0, 
the field equations (8.1) imply the local conservation the total energy and 
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momentum, 

V^T'^^ = . (8.2) 

which can similarly be split into energy and momentum conservation equations 
with respect to a given coordinate system. 

In the case of multiple matter components the total energy-momentum ten- 
sor is the sum of the energy-momentum tensors of the individual fluids, T^'^'^ , 
given in Eq. (4.20). For each fluid the local energy-momentum "conservation" 
equation (4.21), has an energy-momentum transfer 4- vector, Q\^a) right- 
hand-side, which is zero only for non-interacting fluids. However local conser- 
vation of the total energy-momentum imposes the constraint Eq. (4.22). 

We also have at our disposal the equations of motion for specific matter fields, 
such as the Klein-Gordon equation for canonical scalar fields, 9?/, with inter- 
action potential energy U : 

dU , , 

U^i = —. 8.3 

In the following we equate terms order by order in a perturbative expansion 
about a homogeneous background spacetime. 



8.1 Background 



The Einstein equations (8.1) give the Friedmann constraint and evolution 
equation for the background FRW universe 

n^^^a^p, (8.4) 
n'^-^a^p + 3P) , (8.5) 

and energy-momentum conservation, Eq. (8.2), gives the continuity equation 

p' ^-3n{p + P) , (8.6) 

where p and P are the total energy density and the total pressure, a prime 
denotes a derivative with respect to conformal time, rj, the scale factor is a, 
and H = a' /a is the conformal Hubble parameter. 

The total density and the total pressure are related to the density and pressure 
of the component fluids by 

X:Pa = P, E^" = P. (8.7) 
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The continuity equation (4.21) for each individual fluid in the background is 
[66] 

= -m (p, + P,) + aQ^ , (8.8) 

where the energy transfer to the a-fluid is given by the time component of the 
energy-momentum transfer vector 

g« = -u^gf,) . (8.9) 

Equation (4.22) imphes that the energy transfer obeys the constraint 

E^« = 0. (8.10) 



Homogeneous scalar flelds in the FRW metric obey the Klein-Gordon equation 

ip'; + 2nipi + a''^^0. (8.11) 
d(pi 

It is sometimes useful to identify the kinetic energy density and (isotropic) 
pressure of each field as 

p, = P, = ia-V/ (8.12) 

The Klein-Gordon equation (8.11) then implies an energy transfer of the form 
given by Eq. (8.8) 

aQi = -ip'^. (8.13) 
where this energy is transferred to the potential energy 

pu = -Pu = U, (8.14) 

and overall energy conservation (8.10) implies 

<5c/ = -E<5/ = a~'C/'- (8.15) 



8.2 First order scalar perturbations 



In the following we discuss the linear constraint and evolution equations for in- 
homogeneous perturbations at first order. We omit the subscript "1" denoting 
the order of the perturbations to avoid unnecessary clutter. 



8.2.1 Einstein equations 

The scalar metric perturbations in an arbitrary gauge are related to matter 
perturbations via the first-order energy and momentum constraints [66,123] 
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m {tP' + Hcf)) - [V' + Ha] = -ATtGa^dp, 

jj;' + n(l)^ -AnGa^ip + P)V 



(8.16) 
(8.17) 



where the total covariant velocity perturbation is given by 

V = v + B, (8.18) 
and V is the total scalar velocity potential (4.5). 

In a specific gauge, such as the spatially flat gauge these can be written in 
terms of the corresponding gauge-invariant quantities. For instance, in the 
spatially flat gauge we have 



S^Vflat - ^VVflat = -47rG'a'5pflat, (8.19) 
n^^^ = -AirGa^p + P)V^,t , (8.20) 

which makes it straightforward to eliminate the metric variables (p^^^ and 
Cflat = — -Bflat in favour of the energy and momentum in the flat gauge. 

Alternatively in the longitudinal gauge the shear terms are absent and we 
obtain flrst order differential equations for the curvature perturbation 



37i (*' + 7i$) - V^^ = -47rGa25p<,, (8.21) 
^f' + n^^-A7rGa^{p + P)ve. (8.22) 

Typically one then uses these equations to eliminate the density and velocity 
perturbations, Sp^ and V£, in terms of the metric perturbations in the longitu- 
dinal gauge. 

The same energy and momentum constraints can be re-written in terms of 
gauge invariant variables to give expressions for the curvature perturbation 
in the uniform-density gauge (7.61) and the comoving curvature perturbation 
(7.46), respectively, in terms of the longitudinal gauge metric perturbations 
(7.1) and (7.2): 



nil nj2 ul u2 



n 



n 



(8.23) 
(8.24) 



These can be combined to give the gauge-invariant generalisation of the New- 
tonian Poisson equation 

V'* = -3 (n' - Tf) (C + ^) = 47rG'a'5pcom , (8.25) 
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relating the longitTidinal gauge curvature perturbation (7.2) to the comoving 
density perturbation (7.48). We see that the comoving density perturbation is 
suppressed relative to the metric perturbation ^' on large scales, and that ( 
and —TZ coincide in the large scale hmit so long as ^ is finite in this hmit ^ . 

The perturbed Einstein equations at first order also yield two evolution equa- 
tions for the scalar metric perturbations 

ip" + 2nip' + net)' + (^n' + n'')(f)^ AirCa'' (^sp + ^ v'n^ , (8.27) 

a + 2na + ^ - = STiGa'^U , (8.28) 

where IT is the scalar part of the (tracefree) anisotropic stress, defined in 
Eq. (4.10). 

Equation (8.28) in a general gauge can be interpreted as the evolution equation 
for the scalar shear, but in the longitudinal gauge it becomes a constraint 
equation for the gauge-invariant perturbations $ and ^, defined in Eqs. (7.1) 
and (7.2), 

^ - <^ = STiGa^U , (8.29) 
and hence we have ^' = $ in the absence of anisotropic stresses. 

Equation (8.27) then provides a second-order evohition equation for the metric 
perturbation in the longitudinal gauge driven by isotropic pressure: 

+ 37^*' + (2n' + 7^2) * = AnGa^SP . (8.30) 

For adiabatic perturbations we can relate the pressure to the density, SP = 
c^Sp where cj. is the adiabatic sound speed, in which case (8.30) and (8.16) 
yield a closed second-order differential equation [123] 

+ 3(1 + c^)7i*' + [2n' + (1 + 3cl)7f - c^V^]* = . (8.31) 
8.2.2 Energy and momentum conservation 

Energy-momentum conservation gives evolution equations for the perturbed 
energy and momentum 

^ Note that Eq. (8.24) shows that the variable denoted ( in the review by 
Mukhanov, Feldman and Brandenberger: 

Cmfb ^ $ - \^,_!fl2 > (8-26) 
coincides with the comoving curvature perturbation when * = 
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V + 37i(5p + 5P)-3(p + P)V'' + (p + P)V'(l^ + a) = 0, (8.32) 

v' + {i- 3(Pjnv + + -^-^ (^sp + ^ v'n) = o , (8.33) 



p 

where Cg is the adiabatic speed of sound, defined as 



c!^j. (8.34) 



Prom the momentum conservation equation in the total matter gauge, such 
that V — 0, we see that the acceleration is proportional to the pressure per- 
turbation: {p + P)(^ = SP + (2/3) V^n. Alternatively, for pressureless, non- 
interacting dust we have {aV)' + a(f) = and hence the scalar velocity potential 
redshifts as oc 1/a in a synchronous gauge. 

Re-writing the energy conservation equation (8.32) in terms of the curvature 
perturbation on uniform-density hypersurfaces, ^ in (7.61), we obtain the im- 
portant result 

C = - , (8.35) 

p + P 

where (5Pnad is the non-adiabatic pressure perturbation, defined in (9.5), and 
S describes the divergence of the velocity in the longitudinal gauge, Eq. (7.8) 
or, equivalently, the scalar shear along comoving worldlines [96] 

Sy = lv2(V + a) = ^VV,. (8.36) 

Thus the curvature perturbation on uniform-density hypersurfaces is constant 

for adiabatic perturbations on large scales when the shear of comoving world- 
lines becomes negligible. This follows directly from local energy conservation 
and holds independently of the gravitational field equations [176,159,22,33]. 

Using the definition of ^ and ^ in Eqs. (7.61) and (7.2) and the constraint 
equation for the comoving density perturbation (8.25), we have 

^=^(C + *) + .7.=h^l^l (8.37) 

Thus we see that ( is constant for adiabatic perturbations (5Pnad = 0) on 
super-Hubble scales {k/7i -C 1), so long as ^' remains finite. This makes 
C a convenient variable to characterise the primordial density perturbation 
on super-Hubble scales, either during a period of inflation in the very early 
universe, or in the subsequent radiation dominated era. This is an excellent 
approximation throughout reheating at the end of inflation and the subsequent 
radiation era on scales relevant for observations of degree-scale anisotropics in 
the cosmic microwave background and large-scale galaxy surveys [84]. 
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Conversely, local variations in the pressure leading to a non-adiabatic pressure 
perturbation, (5Pj,ad) will lead to a change in the curvature perturbation ^ on 
super- Hubble scales [116,50,89]. This mechanism is at the heart of the curvaton 
scenario for the origin of large scale structure in the universe [48,95,118]. 



8.2.3 Multiple fluids 

The perturbations in the total energy-momentum can be related to the per- 
turbations of individual fluids by 

Y^Spa^Sp, Y.^P-^5P, Ena = n, (8.38) 

a Q. Q. 

and 

^ = E^^^7, (8.39) 

where 5pa and 5Pa are the perturbed energy density and the perturbed pres- 
sure of the CK-fluid, respectively, and Va is the covariant velocity perturbation 
of the a-fluid defined as 

V^ = v^ + B, (8.40) 
where Va is the scalar velocity potential of the a-fiuid. 

The perturbed energy transfer 4- vector, Eq. (4.21), for individual fluids in- 
cluding terms up to first order, is written as [66,113] 



Q(a)o = -a(5a(l + 0) - aSQa , 

Q(a)i = {fa + aQaV)^, , (8.41) 

and Eq. (4.22) implies that the perturbed energy and momentum transfer 
obey the constraints 

^SQa^O, E/a = 0. (8.42) 

a a 

Note that the momentum transfer, f^, is by convention [66,111] defined with 
respect to the total momentum, V, so is non-zero only if the momentum 
transfer vanishes in the total matter frame {V — 0). 

The perturbed energy conservation equation for a particular fluid, including 
energy transfer, is then obtained by the first-order part of the time-component 
of the perturbed continuity equation (4.21) to give [66,104] 

Va + 3H(5pa + SPa) -3{pa + Pa) ^' + {Pa + Pa)^^ (K. + f^) = oQa^ + dSQa , 

(8.43) 

The momentum conservation equation of the a-fiuid is 
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v + 

+ 



aQc 



Pa + Pc 

1 

pa + Pa 



Va + 



(1 + ci) + (1 - 2,cl)H 



= 0, 



(8.44) 



where = Pa/ Pa IS the adiabatic sound speed of the a-fluid and a^[n„,jj — 
{l/3)6ij'V'^Ila] is the scalar anisotropic stress of that fluid. The total adiabatic 
sound speed, Eq. (8.34), is the weighted sum of the adiabatic sound speeds of 
the individual fluids. 



Pa „2 



.45) 



We recover the evolution equation for the total density perturbation (8.32) 
from Eq. (8.43) by summing over all fluids, using Eq. (8.38) and the constraint 
(8.42). 

Analogous to the curvature perturbation on uniform-total-density hypersur- 
faces, C deflned in Eq. (7.61), we can deflne a gauge-invariant perturbation on 
the uniform- CK-density hyper surf aces 



A _ / nj^Pa 
Ca = --0 - T 
Pa 



(8.46) 



The perturbed energy conservation equation for each fluid can then be written 

as 



Pa 



P'a 



7^. 



/HV aQa ( 5 Pa 5p 
\aj p'a 



(8-47) 

P J Pa V Pa J 



where the intrinsic non-adiabatic pressure perturbation of each fluid is given 

by 

SPintr,a = SPa-clSpa, (8.48) 

the scalar shear along worldlines comoving with the a-fluid is 



-y{<j+Va) , 



(8.49) 



and, extending Eq. (7.46), the curvature perturbation comoving with the a- 
fluid is 

7^« = V + n{va + B) . (8.50) 
Thus we see that C^a is constant on large scales for adiabatic perturbations of 
any perfect fluid, with 5Pintr,a = 0, whose energy is conserved, Qa = [176]. 
In fact we shall show later that is constant even in the presence of energy 
transfer, Qa 7^ 0, so long as that energy transfer is adiabatic [111]. 
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8.2.4 Multiple fields 

If we consider scalar fields with Lagrangian density 

1 TV 

C = -Ui^,, . . . , <^^) - - ^ 5^V/,M<^v > (8-51) 

^ 1=1 

and minimal coupling to gravity, then the total energy, pressure and momen- 
tum perturbations are given by 

6p = J2 [«" {S^'i - + UiSifi] , (8.52) 

5P^Y1 [«~ V/ {^¥^'1 - <pM - Ui5<pj] , (8.53) 
I 

{p + P){v + B),i = -Y^ a-'ip'j6ipi,i , (8.54) 

where Uj = dU/dipi. These then give the gauge-invariant comoving density 
perturbation (7.48) 



Spm = ^ [if', {5ip'j - ipM - {ip'; - Wi)Sipi] . (8.55) 
I 

The comoving density is sometimes used to represent the total matter per- 
turbation, but for a single scalar field it is proportional to the non-adiabatic 
pressure (9.5): 

5Pnad = -^^5pcom. (8.56) 

From the Einstein constraint equation (8.25) this will vanish on large scales 
{k/aH ^ 0) if \[' remains finite, and hence single scalar field perturbations 
become adiabatic in this large-scale limit, even without assuming slow- roll 
[52,104]. 

The anisotropic stress, Tiij, vanishes to linear order for any number of scalar 
fields minimally coupled to gravity. 

The first-order scalar field perturbations obey the wave equation 



J 

(8. 



where the terms on the right-hand-side represent the effect of metric pertur- 
bations at first-order (sometimes called the gravitational back-reaction). 
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The scalar field wave equation in a perturbed FRW cosmology can most easily 
be written in closed form in terms of the field perturbations in the spatially 
fiat gauge, defined in Eq. (7.23), which in the multi- field case have the gauge 
invariant definitions 

5v?flat/ = ^^^/ + f^^- (8.58) 

Note that since the scalar field can be thought of as a potential for the 
4-velocity, this variable is a rescaling of the curvature perturbation on the 
comoving-orthogonal hypersurface for each field 

<^^flat7 = ■ (8-59) 

Using the Einstein equations to eliminate the metric perturbations on the 
right-hand-side of Eq. (8.57) yields [146,64,128,164] 



(8.60) 

The effect of gravitational coupling is now evident due to the terms propor- 
tional to Newton's gravitational constant. It is also evident that this gravita- 
tional coupling vanishes at first order for fields whose time derivative vanishes 
in the background solution, which is why at lowest order in a slow-roll approx- 
imation during inflation one can neglect the gravitational coupling. 

In the next section we will discuss how the coupled equations for multiple fields 
may be partially decoupled by identifying the adiabatic and isocurvature field 
perturbations on large scales. 



8.3 First order vector perturbations 



The divergence-free part of the 3-momentum [see Eqs. (2.11), (4.5) and (4.13)] 

^g^ = (P + ^)(^vec^-^^), (8.61) 

obeys the momentum conservation equation 

5(i, + AWqi = -V^Ui , (8.62) 

where the vector part of the anisotropic stress, Eq. (4.10), is given by 0^9(^11^). 
The gauge-invariant vector metric perturbation is then directly related to the 
divergence-free part of the momentum via the constraint equation 

{f; + Si) = -167rGa^5qi . (8.63) 
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Thus the Einstein equations constrain the gauge-invariant vector metric per- 
turbation, FI + Si, to vanish in the presence of only scalar fields, for which the 
divergence-free momentum necessarily vanishes. 

Equation (8.63) shows that vector metric perturbations can be supported only 
by divergence-free momenta, but even then equation (8.62) shows that the 
vector perturbations are redshifted away by the Hubble expansion on large 
scales unless they are driven by an anisotropic stress. 

8.4 First order tensor perturbations 

There is no constraint equation for the tensor perturbations as these are the 
free gravitational degrees of freedom (gravitational waves). The spatial part 
of the Einstein equations yields a wave equation 

h'lj + 2nhlj - V^hij = STrGa^n^ , (8.64) 
where H^J^-* is the transverse and tracefree part of the anisotropic stress (4.10). 

We can decompose arbitrary tensor perturbations into eigenmodes of the spa- 
tial Laplacian, V^Cj^ = —{k'^/a'^)eij, with comoving wavenumber k, and scalar 
amplitude h{t): 

= hit)e'i^'''\x) , (8.65) 
with two possible polarisation states, + and x. In the absence of any such 
anisotropic stress, e.g., in the presence of scalar fields and perfect fiuids, the 
amplitude, defined in Eq. (8.65), of the tensor metric perturbation with co- 
moving wavenumber, k, obeys the wave equation for a massless scalar field 
(8.57) in an unperturbed FRW metric. 

h" + 2nh' + k^h^O. (8.66) 



9 Adiabatic and entropy perturbations 

9. 1 Multiple fluids 

We will refer to primordial perturbations as the perturbations at the epoch 
of primordial nucleosynthesis. The abundances of the light elements provide 

constraints on the matter content and expansion rate of the universe at this 
epoch, so we will assume that the material content of the universe is known 
(photons, neutrinos, baryonic matter and cold dark matter) and the gravita- 
tional laws are described by general relativity. This is expected to be some 
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time after an early inflationary epoch when the perturbations originated as 
vacuum fluctuations on much smaller scales. By the time of primordial nucle- 
osynthesis, the scales responsible for the large scale structure of our observable 
universe today were far outside the Hubble scale and well described by large 
scale limit. 

In the standard Hot Big Bang the entropy density of the universe is domi- 
nated by the number of relativistic particles and for most of the history of 
the universe it is proportional to the number of photons, s = l.SgsTi^, where 
Qs is the effective number of light species [70]. In particular the perturbed 
baryon-entropy ratio ub/s (assuming gs remains constant) is given by [116] 

_S{nB/n^) Sub 5n^ 

^B = -, — • 

nB/n^ ub 

Written in terms of the energy density of photons, oc n^^^, and baryons, 
Pb oc Ub, Bi the time of primordial nucleosynthesis this becomes 

g^^SpB_35p, 
Pb 4: p^ 



More generally, density perturbations in an n-componcnt system can be de- 
composed into an overall density perturbation and n — 1 relative density per- 
turbations between the different components. The overall density perturbation 
is naturally gauge-dependent, however the gauge transformation rule for the 
linear density perturbation (6.18) suggests a natural gauge-invariant deflnition 
of the relative density perturbation at flrst-order, 

5pij (X 5pi - ^5pj , (9.3) 
Pj 

corresponding to the density perturbation of fluid / on surfaces of uniform den- 
sity of the fluid J. Comparing this expression with the conventional deflnition 
of the primordial baryon-entropy perturbation (9.2) suggests a gauge-invariant 
deflnition of the relative perturbation between any two fluids [113,111] 

5„^3W(^-^)^ 3(0-0). (9.4) 
\ Pj Pi ) 

which reduces to baryon-entropy perturbation (9.2) for Sb = Sb'^- Hence we 
refer to Su as the relative entropy perturbation for two fluids. It is the correct 
generalisation of the the entropy perturbation deflned in Ref. [66] to the case 
of non-interacting fluids. 
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The non-adiabatic pressure perturbation is given by 

SP^^ = SP-^Sp. (9.5) 
P 

For a detailed recent discussion of the non-adiabatic pressure see Ref. [37]. 

In a system of more than one fluids the total non-adiabatic pressure pertur- 
bation, 5Pnad) ms-y be further split into two parts [66], 

5P„ad = SPintr + SP^el • (9.6) 

The first part is due to the intrinsic entropy perturbation of each fluid 

<^Pintr = E<^^nad,a, (9.7) 
a 

where the intrinsic non-adiabatic pressure perturbation of each fluid was given 
in Eq. (8.48). The second part of the non-adiabatic pressure perturbation (9.6) 
is due to the relative entropy perturbation S^p between different fluids (9.4) 

<^Prel = E P'Jp (4 - 4) Sap • (9.8) 



In analogy with the non-adiabatic pressure perturbation for each fluid (8.48), 
we can identify an intrinsic non-adiabatic part of the energy transfer pertur- 
bation [113,111] that appears in the perturbed energy conservation equation 
for each fluid (8.43) 

SQintr,a = SQa - ■ (9.9) 

Pa 

This is automatically zero if the local energy transfer is a function of 
the local density Pa so that 6Qa = {dQa/dpa)Spa, just as the intrinsic non- 
adiabatic pressure perturbation (8.48) vanishes when SPa = {dPa/dpa)5pa- 

We can also identify in Eq. (8.43) a relative non-adiabatic energy transfer 



An _^ 'H'-n' (5p^ 5p' 

0^rel,a ^a ^ 

Qc 



,^^^P',Sap, (9.10) 

due to the presence of relative entropy perturbations whenever the background 
energy transfer is non-zero, Qa 7^ 0. 

The perturbed energy conservation Eq. (8.47) for each fluid can then be writ- 
ten as 
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Pa Pa •^'^ 



Pa 



(9.11) 



where the non-adiabatic pressure perturbation of each fluid is given by Eq. (8.48) 
and the non-adiabatic energy transfer is given by Eqs. (9.9) and (9.10). We 
thus see that C,ct is constant for adiabatic perturbations in the large scale limit 
where the shear of comoving worldlines, defined in Eq. (8.49), vanishes 
[111]. 

More generally [96] one finds a conserved perturbation whenever there is a 
local conservation equation of the form dy/dr — —9f{y), where 9 is the local 
expansion rate and r the proper time along comoving worldlines. When one 
integrates this conservation equation one finds the local logarithmic expansion 
as a function of y: 

Thus the difference 5F = F^ys) — F{yA) evaluated along different worldlines 
remains a fixed constant of integration if one evaluates 5F on constant time 
hypersurfaces separated by a fixed expansion, and spatially fiat hypersurfaces 
provide a suitable time-slicing on large scales. The classic example of such a 
conserved quantity is the curvature perturbation on uniform density hypersur- 
faces, C defined in Eq. (7.61), which is conserved on large scales when P — P{p) 
and thus the perturbations are adiabatic. But as we have seen the same result 
holds for the curvature perturbation on uniform- density hypersurfaces for any 
fluid whose pressure perturbation and energy transfer is adiabatic. 

Thus the matter isocurvature perturbation (9.2) is constant on large scales if 
there is negligible energy transfer between non-relativistic matter and radia- 
tion. However at even higher energies wc can still define a conserved pertur- 
bation associated with conserved baryon number density [97] 



5nB 
riB 



-ip (9.13) 



so long as we have a conserved quantum number associated with baryon num- 
ber. The observed stability of the proton implies that baryon number is con- 
served up to very high energy, possibly the GUT scale, and thus it should be 
possible to relate any primordial baryon isocurvature perturbation to physics 
at very high energies. 
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9.2 Multiple fields 



In the background FRW cosmology driven by multiple scalar fields it is possible 
to identify an adiabatic direction along the background trajectory in field space 

ri^4^. (9.14) 



The background solution, even in the presence of multiple fields, can then 
be described in terms of an effective single field, r, obeying the usual Klein- 
Gordon equation 

r" + 2nr' + Ur = 0. (9.15) 

where r' = J^i'f'if'i f-nd Ur = Yli^iUi- However field perturbations need 
not follow this background trajectory and we encounter quaUtatively different 
behaviour from that in the single field case when we consider inhomogeneous 
perturbations about the background trajectory. 

In analogy with our treatment of fluid perturbations one can identify adiabatic 
and isocurvature field perturbations in a cosmology with more than one scalar 
field. Indeed even a single scalar field, </?, can support a non-adiabatic pressure 
perturbation, given in Eq. (8.56). We refer to this as the intrinsic pressure 
perturbation for the field. However this intrinsic non-adiabatic perturbation 
for a single field is proportional to the comoving density perturbation and 
thus vanishes in a scalar field dominated universe on large scales according to 
Eq. (8.25), leaving effectively only adiabatic perturbations in this large scale 
limit. 

For multiple fields we can perform a local rotation in field space to identify 
the adiabatic part of arbitrary perturbations along the background trajectory 
[109,52] 

5r = Y,fi^'Pi- (9-16) 

The generalisation to non-canonical fields with arbitrary metric in field space 
is given in Refs. [55,41]. 

Field perturbations orthogonal to the adiabatic field are isocurvature field 
perturbations, or entropy perturbations in analogy with the fiuid density per- 
turbations (9.2), 

dsij(x— — . (9.17) 

Note that the adiabatic field perturbation is naturally gauge- dependent, whereas 
the isocurvature field perturbations are gauge-independent at first order. For 
simplicity we will consider the case of two scalar fields where the direction in 
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^1 



Fig. 1. r and s are used to denote the instantaneous adiabatic and entropy fields 
respectively, along and orthogonal to the curved background solution in field space. 

field space is given by the angle 6*, see Figure 1, so that f = (cos 6*, sin^). We 
then have 



5r = cos 96ipi + sin 96ip2 , 
6s = — sin 9S(pi + cos 96ip2 



(9.18) 
(9.19) 



and we will work with the adiabatic field perturbation in the spatially fiat 
gauge 



ri 



(9.20) 



Note that 5rflat is thus proportional to the total comoving curvature pertur- 
bation 

n = ^+'^6r = ^6r^^, . (9.21) 
since the adiabatic field, r, is the potential for the total velocity, oc dr /dx^. 

The adiabatic and isocurvature field perturbations obey the coupled evolution 
equations [52] 



flat 



A;2 + a^U„ - 9''- 



2{9'6s) 



SttG /aV' 



n 



at 



+ 



n 

9' 



9' 5s. 



6s" + 2ms' +{k'' + a^Uss + 3^'2) 5s = 



where 



(9.22) 
(9.23) 



Vrr = (cos' 9)Uu + (siu 29)Uu + (sin' 9)U22, 
Uss = (sin' 9)Uu - (sin 29)Ui2 + (cos' 9)U22. 



(9.24) 
(9.25) 
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We can identify a purely adiabatic mode where 5s = and remains zero 
on large scales. However a non-zero isocurvature perturbation appears as a 
source term in the perturbed inflaton equation (9.22) whenever 6' ^ and the 
inflaton trajectory is curved in field space. Note that 9' is given by [52] 

e'^-'L^^ (9.26) 

where Ug — (cos 9)U2 — (sin 9)Ux is the potential gradient orthogonal to the in- 
flaton trajectory in fleld space. In the slow-roll approximation the background 
fleld always follows the potential gradient so the adiabatic-isocurvature cou- 
pling is suppressed in this slow-roll limit the integrated effect of isocurvature 
field perturbations on the adiabatic field perturbation cannot in general be 
neglected. 

Equation (9.22) shows that the isocurvature perturbation 5s works as a source 
term for the adiabatic curvature perturbation. This is in fact clearly seen if 
we take the time derivative of the comoving curvature perturbation (9.21): 

n' = :^^^k-^+'^9'5s. (9.27) 

Therefore TZ (or equivalently C) is not conserved even in the large-scale limit 
in the presence of an isocurvature field perturbation, 5s, with a non-straight 
trajectory in field space {9' ^ 0). 

By contrast, the solution for the isocurvature field perturbation is independent 
of any initial adiabatic perturbation on large-scales. The adiabatic perturba- 
tion provides a source term for the isocurvature field only through the spatial 
gradient of the longitudinal gauge metric potential, ^, which rapidly becomes 
negligible on super-Hubble scales during slow-roll inflation. 



10 Perturbations from inflation 



The standard hot big bang model of cosmology has a major shortcoming in 
that there is no causal explanation for the existence of primordial density per- 
turbations on super-Hubble scales during the radiation dominated era. The 
CMB, and in particular the acoustic peaks in the temperature and polari- 
sation anisotropics seen by the WMAP sateUite [74] provide strong evidence 
that these primordial density perturbations do exist on scales much larger 
than the causal horizon at early times. The detailed distribution of primordial 
inhomogeneities is left as an unexplained initial condition in the standard hot 
big bang. 
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The primary success of inflation [155,56] is to give a model for the origin of the 
primordial density perturbations from vacuum fluctuations during a period 
of accelerated expansion at very early times. This relies on speculative and 
uncertain physics - in particular it requires some form of energy density with 
negative pressure in general relativity - but the unexpected discovery that the 
universe is accelerating today appears to show that cosmological inflation does 
happen. Zero-point vacuum fluctuations of any light, weakly-coupled scalar 
field will be stretched up to super-Hubble scales during infiation and leave an 
approximately scale-invariant and Gaussian distribution of perturbations on 
large scales [122,58,156,57,9,119,145,120]. 

The simplest model for inflation is that it is driven by a vacuum (potential) 
energy density which is a function of one or more scalar fields. We can describe 
the homogeneous FRW solution using the inflaton field, r, which describes the 
evolution along the trajectory, r defined in Eq. (9.14), in a possibly multi- 
dimensional field space. For a sufficiently fiat potential the evolution can be 
well- described by the slow-roll approximation which assumes that the energy 
density is potential dominated 

^ ^^U{r) . (10.1) 
This is equivalent to requiring that the first slow-roll parameter is small: 



<1. (10.2) 



Note that the condition for accelerated expansion requires e < 1. In terms of 
the potential we have 

We also assume the universe is overdamped, such that 

SHr' ~ -a^Ur . (10.4) 

This implies that we can neglect the decaying mode of the overdamped system 
and we have a unique trajectory in field space for a single field. This is a self- 
consistent approximation when the infiaton field is light compared with the 
Hubble scale, which requires that the second slow-roll parameter is small: 

\r)rr\ < 1 where rjrr = ■ (10.5) 



In the two-field model described in section 9.2, allowing for non-adiabatic 
perturbations, we can define three parameters [175], r]rr, rjrs and rjss, describing 
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the curvature of the potential, where in general we have 



The background slow-roll solution is described in terms of the slow-roll pa- 
rameters by 

r^~^eC/, H-^^-qrs, (10.7) 
while the perturbations on large scales (comoving wavenumber k <^7i) obey 



H ^^rflat ^ (2e - r}rr) Sr^ax - 2r)rsSs , 

H-^Ssc^-r]ssSs, (10.8) 

where we neglect spatial gradients. Although C/g ~ at lowest order in slow- 
roll, this does not mean that the inflaton and entropy perturbations decouple. 

6 given by Eq. (10.7) is in general non-zero at first-order in slow-roll and 
large-scale entropy perturbations do affect the evolution of the adiabatic per- 
turbations when r]„s 7^ 0. 

While the general solution to the two second-order perturbation equations 
(9.22) and (9.23) has four independent modes, the two first-order slow-roll 
equations (10.8) give the approximate form of the squeezed state on large 
scales. This has only two modes which we can describe in terms of the dimen- 
sionless comoving curvature and isocurvature perturbations: 

7^=^5rflat, S='^5s. (10.9) 

The normalisation of TZ coincides with the standard definition of the comoving 
curvature perturbation, Eq. (7.46). The normalisation of the dimensionless 
entropy during inflation, S, is chosen here coincide with Ref. [175]. It can be 
related to the non-adiabatic pressure perturbation (9.5) on large scales 

SPn.d--eVrs^S. (10.10) 



The slow-roll approximation can provide a useful approximation to the instan- 
taneous evolution of the flelds and their perturbations on large scales during 
slow-roll inflation, but is not expected to remain accurate when integrated 
over many Hubble times, where inaccuracies can accumulate. In single-field 
infiation the constancy of the comoving curvature perturbation after Hubble 
exit, which does not rely on the slow- roll approximation, is crucial in order to 
make accurate predictions of the primordial perturbations using the slow-roll 
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approximation only around Hubble crossing. In a two-field model we must 
describe the evolution after Hubble exit in terms of a general transfer matrix: 





(10.11) 



On large scales the comoving curvature perturbation still remains constant 
for the purely adiabatic mode, corresponding to 5 = 0, and adiabatic per- 
turbations remain adiabatic. These general results are enough to fix two of 
the coefficients in the transfer matrix, but Tus and T55 remain to be deter- 
mined either within a given theoretical model, or from observations, or ideally 
by both. The scale-dependence of the transfer functions depends upon the 
infiaton-entropy coupling at Hubble exit during inflation and can be given in 
terms of the slow-roll parameters as [175] 



d 

q^^Ttzs = '2r]rs + (2e - r]rr + Vss)Tns , 

^^Tss = {2e-7]rr + Vss)Tss. (10.12) 



10.1 Initial power spectra 



The expectation value of the fluctuations of a homogeneous field are given by 

{5ipiik,)5ipiik2)) = (27r)^5^(ki + k2)Ps^,{h) , (10.13) 

where angle brackets denote the ensemble average (i.e., the average over in- 
finitely many realisations of the field, which is equivalent to taking the spatial 
average in an infinite space). The dimensionless power spectrum of the field 
(equivalently the variance of the field per logarithmic range of k) is given by 

T'^^ii^) = T^^P^^i^^) ■ (10-14) 



For weakly-coupled, light fields we can neglect interactions on wavelengths 
below the Hubble scale, so that vacuum fiuctuations correspond to [86] 

^5^.W = ^- (10.15) 
This gives rise to the power spectrum of field fiuctuations on the Hubble scale 
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{k — H — aH) during inflation given by 



v^, . (I); , (10.16) 

where we use a * to denote quantities evaluated at Hubble-exit. If a field 
has a mass comparable to the Hubble scale or larger then the vacuum fluc- 
tuations on wavelengths greater than the effective Compton wavelength are 
suppressed. In addition fluctuations in strongly interacting fields may develop 
correlations before Hubble exit. But during slow- roll inflation the correlation 
between vacuum fluctuations in weakly coupled, light fields at Hubble-exit is 
suppressed by slow-roll parameters. This remains true under a local rotation 
in fields space to another orthogonal basis such as the instantaneous inflaton 
and entropy directions (9.18) and (9.19) in field space. 

The curvature and isocurvature power spectra at Hubble-exit are given by 

ao.ir, 

while the cross-correlation is first-order in slow- roll [55,30], 

Cnslc^-2C7^rsVnl , (10.18) 

where C — 2 — la2 — ^ ^ 0.73 and 7 is the Euler number. The normalisation 
chosen for the dimensionless entropy perturbation in Eq. (10.9) ensures that 
the curvature and isocurvature fluctuations have the same power at horizon 
exit [175]. The spectral tilts at horizon-exit are also the same and are given 

by 

AnnU ^ ^nsU ^ -6e + 2r]rr . (10.19) 
where Anx = dlnVx/dlnk. 

The tensor perturbations (8.64) are decoupled from scalar metric perturba- 
tions at first-order and hence the power spectrum has the same form as in 
single field inflation. Thus the power spectrum of gravitational waves on super- 
Hubble scales during inflation is given by 



16H^ 128 

-p„ r ~ 1_ do 20) 



and the spectral tilt is 

AnT\*^-2e. (10.21) 
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10.2 Primordial power spectra 



The resulting primordial power spectra on large scales can be obtained simply 
by applying the general transfer matrix (10.11) to the initial scalar perturba- 
tions. The scalar power spectra probed by astronomical observations are thus 
given by [175] 



Vn^{l+Tls)Vn\* (10.22) 
Vs = TlsVnU (10.23) 
Cns^T-^TssVnU. (10.24) 

Note that the primordial curvature and isocurvature perturbations from infla- 
tion are in general correlated [76] (see also [28,168,4]). The cross-correlation 
can be given in terms of a dimensionless correlation angle: 



We see that if wc can determine the dimensionless correlation angle, 0, from 
observations, then this determines the off-diagonal term in the transfer matrix 

T7j5 = cote, (10.26) 

and we can in effect measure the contribution of the entropy perturbation 
during two-field inflation to the resultant curvature primordial perturbation. 
In particular this allows us in principle to deduce from observations the power 
spectrum of the curvature perturbation at Hubble-exit during two-field slow- 
roll inflation [175]: 

Vn\*=Vn&^^^Q- (10.27) 

The scale-dependence of the resulting scalar power spectra depends both upon 
the scale-dependence of the initial power spectra and of the transfer coeffi- 
cients. The spectral tilts are given from Eqs. (10.22-10.24) by 



/^nn = AriTel* + H:^{dTns/dU) sin 20 , 

An5 = AriTel* + 2H;\d\nTss I dU) , (10.28) 
Anc = Ar^7^|. + H;^ [{dT-us/dU)ts.nQ + {d\nTss/dU)] , 

where we have used Eq. (10.26) to eliminate T-jis in favour of the observable 
correlation angle G. Substituting Eq. (10.19) for the tilt at Hubble-exit, and 
Eqs. (10.12) for the scale-dependence of the transfer functions, we obtain [175] 
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Ann^-{G - 4cos^ e)e 

+2 (rjrr sin^ © + 2rirs sin © cos © + rjss cos^ ©^ , 

Ans--2e + 2r]ss, (10.29) 
Arte ^ — 2e + 2r]ss + 2r]rs tan © . 

Although the overall amplitude of the transfer functions are dependent upon 
the evolution after Hubble-exit and through reheating into the radiation era, 
the spectral tilts can be expressed solely in terms of the slow-roll parameters 
at Hubble-exit during inflation and the correlation angle, 0, which can in 
principle be observed. 

If the primordial curvature perturbation results solely from the adiabatic in- 
fiaton field fiuctuations during infiation then we have T-j^s = in Eq. (10.22) 
and hence cos© = in Eqs. (10.29), which yields the standard single field 
result 

Ann - -6e + 2rirr ■ (10.30) 

Any residual isocurvature perturbations must be uncorrelated with the adia- 
batic curvature perturbation (at first-order in slow-roll) with spectral index 

Ans ~ -2e + 2riss ■ (10.31) 

On the other hand, if the observed primordial curvature perturbation is pro- 
duced due to some entropy field fiuctuations during infiation, we have T-jis 3> 1 
and sin © ~ 0. In a two-field infiation model any residual primordial isocurva- 
ture perturbations will then be completely correlated (or anti-correlated) with 
the primordial curvature perturbation and we have 

AriTz ~ Anc ~ Ans ^ -2e + 2r)ss ■ (10.32) 

The gravitational wave power spectrum is frozen-in on large scales, indepen- 
dent of the scalar perturbations, and hence 

Vt^VtU- (10.33) 

Thus we can derive a modified consistency relation [86] between observables 
applicable in the case of two-field slow-roll infiation: 

^ ~ -SAriT sin^ © . (10.34) 

This relation was first obtained in Ref. [10] at the end of two-field infiation, 
and verified in Ref. [169] for slow-roll models. But it was realised in Ref. [175] 
that this relation also applies to the primordial perturbation spectra in the 
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radiation era long after two-field slow-roll inflation has ended and hence may 
be tested observationally. 

More generally, if there is any additional source of the scalar curvature per- 
turbation, such as additional scalar fields during inflation, then this could give 
an additional contribution to the primordial scalar curvature spectrum with- 
out affecting the gravitational waves, and hence the more general result is the 
inequality [158,146]: 

^ < -SAnx sin^ 9 . (10.35) 

This leads to a fundamental difference when interpreting the observational 
constraints on the amplitude of primordial tensor perturbations in multiple 
inflation models. In single field inflation, observations directly constrain r = 
[Pt/VuI* and hence, from Eqs. (10.17) and (10.20), the slow-roll parameter e. 
However in multiple field inflation, non-adiabatic perturbations can enhance 
the power of scalar perturbations after Hubble exit and hence observational 
constraints on the amplitude of primordial tensor perturbations do not directly 
constrain the slow-roll parameter e. 

Current CMB data alone require r < 0.55 (assuming power-law primordial 
spectra) [74] which in single-field models is interpreted as requiring e < 0.04. 
But in multiple field models e could be larger if the primordial density per- 
turbation comes from non-adiabatic perturbations during inflation. 



11 Non-line£ir evolution and non-Gaussianity 



A powerful technique to calculate the primordial curvature perturbation re- 
sulting from many inflation models, including multi-field models, is to note 
that the curvature perturbation C, defined in Eq. (7.61) can be interpreted as 
a perturbation in the local expansion [156,157,146,98] 

C = 5N, (11.1) 

where 5N is the perturbed expansion to uniform-density hypersurfaces with 
respect to spatially flat hypersurfaces, which is given to first-order by 

C = -//^, (11.2) 

where 5pflat must be evaluated on spatially flat (■^ = 0) hypersurfaces (see 
Sections 3.2 and 7.6). 

An important simplification arises on large scales where anisotropy and spa- 
tial gradients can be neglected, and the local density, expansion, etc., obeys the 
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same evolution equations as in a homogeneous FRW universe [144,146,147,176,96,136,98]. 
Thus we can use the homogeneous FRW sohitions to describe the local evolu- 
tion, which is known as the "separate universe" approach [144,146,147,176,136]. 
In particular we can evaluate the perturbed expansion in difTerent parts of the 
universe resulting from different initial values for the fields during inflation 
using the homogeneous background solutions [146] . 

In the slow-roll approximation the integrated expansion on super- Hubble scales 
from some initial spatially flat hypersurface up to a subsequent fixed density 
hypersurface (say at the epoch of primordial nucleosynthesis) is some function 
of the local field values on the initial hypersurface, N[(pi\^). More generally 
we expect this to hold whenever we can neglect the decaying mode for the field 
perturbations on super-Hubble scales. The resulting primordial curvature per- 
turbation on the uniform-density hypersurface is then 

C = E^^Wt, (11-3) 
I 

where Ni = dN/dtpi and 5<^/flat is the field perturbation on some initial 
spatially-flat hypcrsurfaces during infiation (8.58). In particular the power 
spectrum for the primordial density perturbation in a multi-field infiation can 
be written (at leading order) in terms of the field perturbations after Hubble- 
exit as 

I 

This approach is readily extended to estimate the non-linear effect of field per- 
turbations on the metric perturbations [147,96,94,98]. We can take Eq. (11.1) 
as our definition of the non-linear primordial curvature perturbation, C,, so 
that in the radiation dominated era the non-linear extension of Eq. (11.2) is 
given by [98] 

where pflat(^,x) is the perturbed (inhomogeneous) density evaluated on a spa- 
tially fiat hypersurface and po(^) is the background (homogeneous) density. See 
Refs. [136,138,139,78-80] for alternative approaches to define the non-linear 
extension of C,. 

This non-linear curvature perturbation as a function of the initial field fiuctu- 
ations can simply be expanded as a Taylor expansion [94,150,31,151] 

I ^ I,J " I,J,K 

(11.6) 

where we now identify (11-3) as the leading-order term. 
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We expect the field perturbations at Hubble-exit to be close to Gaussian for 
weakly coupled scalar fields during inflation [103,149,137,150,152]. In this case 
the bispectrum of the primordial curvature perturbation at leading (fourth) 
order, can be written as [49,73] 



Sc(ki, k2, kg) = -fNL [Pdkl)Pdk2) + Pc{k2)Pc{h) + Pc{k-i)Pc{ki)]{n.7) 

where P^ik) = 27r^P^(A;)/A;^, and the dimensionless non-linearity parameter is 
given, using the 5A^-formalism, by^ [94] 

- 6 (NcNCf ■ ^^'-^^ 

Similarly to the bispectrum, the connected part of the trispectrum in this case 
can be written as [151,31] 



rc(ki,k2,k3,k4)=r^L [Pd\ki + ks\)Pdk3)Pc{k4) + (11 perms)] 

54 

+7^9NL [Pcik2)Pcik3)Pcik4) + (3 perms)] . (11.9) 
zo 

where 



NabN^^N^Nc 

2hNABcN''N^N^ 
= 54 {NnN^Y ■ ^^^-^^^ 

The expression for tnl was first given in [2]. Note that we have factored out 
products in the trispectrum with different k dependence in order to define 
the two k independent non-linearity parameters rjv^ and Qnl- This gives the 
possibility that observations may be able to distinguish between the two pa- 
rameters [130]. 



In many cases there is single direction in field-space, x, which is responsible 
for perturbing the local expansion, A'^(x), and hence generating the primordial 
curvature perturbation (11.6). For example this would be the inflaton field 
in single field models of infiation, or it could be a late-decaying scalar field 
[116,89,48] as in the curvaton scenario [95,118]. In this case the curvature 



^ Note, that the factor "5/6" in Eq. (11.8) is a historical convention, due to the 
original definition [73] which was given in terms of the Newtonian potential, which 
on large scales in the matter era is given by $ = —(3/5)^- 
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perturbation (11.6) is given by 

C ^ A^'^Xfiat + In"SxL + In"'SxL + ■■■, (11-12) 

and the non-Gaussianity of the primordial perturbation has the simplest "lo- 
cal" form 

C = Ci + pNLCi + ^9nlC! + ■■■ (11-13) 

where Ci = N'Sxfiax is the leading-order Gaussian curvature perturbation and 
the non-linearity parameters /^vl and qnl, are given by [94,148] 



^5 N" 

6 {N'Y ' 
_25 N'" 
^'^''~54(7V03' 

The primordial bispectrum and trispectrum arc then given by Eqs. (11.7) 
and (11.9), where the non-linearity parameters Jnl and Qnl^ given in Eqs. (11.8) 
and (11.11), reduce to Eqs. (11.14) and (11.15) respectively, and tnl given in 
Eq. (11.10) reduces to 



(11.14) 
(11.15) 



^^^ = Jn^^2^^^l- (11-16) 

Thus T^L is proportional to /^^ (first shown in [130] using the Bardeen poten- 
tial, and in [94] using this notation). However the trispectrum could be large 
even when the bispectrum is small because of the qnl term [130,148]. 

In the case of where the primordial curvature perturbation is generated solely 
by adiabatic fluctuations in the inflaton field, r, the curvature perturbation is 
non-hnearly conserved on large scales [96,98,79] and we can calculate A^', N" , 
N'" , etc, at Hubble-exit. In terms of the slow-roll parameters, we find 

A^' = ?^^ — ^-0(e-^/2) , (11.17) 

N"^-]-Wn,r-2e)^0{l) , (11.18) 

2 mpj e 

N'" ^4^^ Urr - V'rr + Ur) ~ 0{e'/') , (11.19) 



where we have used the reduced Planck mass m|i = {8ttG)~^ and introduced 
the second-order slow-roll parameter = rripiUrUrrr / ■ Hence the non- 
linearity parameters for single field infiation, (11.14) and (11.15), are given 

by 
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fNL = 7:{Vrr-2e), (11.20) 
D 

9NL = ^ {2erjrr - 2r^l + C) ■ (11-21) 

with tnl given by Eq. (11.16). Although there are additional contributions to 
the primordial bispectrum and trispectrum coming from the intrinsic non- 
Gaussianity of the field perturbations at Hubble-exit, these are also sup- 
pressed by slow-roll parameters in slow-roll inflation. Thus the primordial 
non-Gaussianity is likely to be too small to ever be observed in the conven- 
tional infiaton scenario of single-field slow-roll inflation [1,103]. Indeed any 
detection of primordial non-Gaussianity /atl > 1 would appear to rule out 
this infiaton scenario. 

However significant non-Gaussianity can be generated due to non-adiabatic 
field fiuctuations. Thus far it has proved difficult to generate detectable non- 
Gaussianity in the curvature perturbation during conventional slow-roll in- 
flation, even in multiple fleld models [140,171], though see Refs. [32,38]. But 
detectable non-Gaussianity might be produced in non-slow-roll [36] or non- 
canonical scalar fleld inflation [3,81], or when the curvature perturbation is 
generated from isocurvature fleld perturbations at the end of inflation [20,93], 
during inhomogeneous reheating [44,69,178,45,72,29], or some time after in- 
flation has ended in the curvaton model [95,99,11,107,148]. 



12 Summciry and outlook 



Linear perturbations have become part of the standard toolbox of modern cos- 
mology. Earlier confusion surrounding apparently different behaviour found in 
different coordinate bases has largely been resolved through the use of vari- 
ables which have gauge-invariant deflnitions. In Section 7 we have emphasised 
the physical interpretation of these gauge-invariant variables. 

The power spectrum of primordial perturbations revealed by the cosmic mi- 
crowave background and large scale galaxy surveys is a powerful probe of 
inflationary models of the very early universe, and a challenge for alternative 
theories. Linear theory enables us to relate the primordial spectra to quantum 
fluctuations in the metric and matter flelds at much higher energies. In the 
simplest single fleld inflation models, it is possible to equate the primordial 
density perturbation with the curvature perturbation, C, deflned in Eq. (7.61), 
during inflation which remains constant on large scales for adiabatic density 
perturbations. More generally, if one allows for non-adiabatic perturbations 
then it becomes necessary to allow for variation in the curvature perturbation, 
even on super-Hubble scales. Nonetheless it is still possible in many cases to 
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identify the primordial density perturbation with a perturbed expansion in 
the SN approach described in Section 11 where the integrated expansion, A^, 
is a function of the local field values on spatially flat slices during inflation. 
As one goes beyond the textbook examples, it becomes necessary to have a 
clear physical definition of the perturbation variables to consistently extend 
the background FRW equations to the inhomogeneous perturbations. 

The new frontier in the study of cosmological perturbations is the study of 
non-linear primordial perturbations, at second-order and beyond. Many of 
the familiar certainties of linear perturbation theory no longer apply. We have 
shown in Section 6 that quantities that were automatically gauge-independent 
at first order (including the non-adiabatic pressure perturbation, anisotropic 
stress, and the tensor metric perturbation) become gauge dependent at second 
order. We have shown in Section 7 that it is possible to use the same method- 
ology to construct gauge-invariant variables at second (and higher) order. A 
variable with an unambiguous physical meaning will have a gauge invariant 
definition. The resulting gauge invariant definitions inevitably become more 
complicated than those at first order and we have only presented explicit def- 
initions at second order for a few cases. Likewise we have not attempted to 
present the second order dynamical equations in a comprehensive manner as 
was done at first order in Section 8. Our aim has been to provide an introduc- 
tion to some of the issues that arise at higher order. Early works on non-linear 
and second order perturbation theory include Refs. [166,131], more recently 
see also Refs. [134,167,129,125-127]. 

Non-linearities allow additional information to be gleaned from the primordial 
perturbations. Much effort is currently being devoted to the study of higher- 
order correlations. Non-Gaussianity in the distribution of primordial density 
perturbations would reveal interactions beyond the linear theory. Such in- 
teractions are minimal (suppressed by slow-roll parameters) in the simplest 
single field inflation model, so any detection of primordial non-Gaussianity 
would cause a major upheaval in our thinking about the very early universe. 
In principle the 6N approach can be easily extended to higher-orders, simply 
by extending the Taylor series for the integrated expansion as a function of the 
fleld perturbations beyond linear order. This enables one to compute higher- 
order correlation functions for ^ as shown in Section 11. The challenge is then 
to develop transfer functions to relate the primordial ( to observables be- 
yond linear order [15,12-14], although large primordial non-Gaussianity (e.g., 
Jnl ^ 1) is expected to dominate over non-linearity in the transfer functions. 

The Klein- Gordon equation in closed form at second order shows that at sec- 
ond order scalar perturbations will also be sourced by terms quadratic in first 
order field perturbations [106,112]. This and second order perturbation theory 
in general, provides an alternative to using the SN formalism in calculating 
the primordial non-gaussianity /jvl [153]. The main advantage of perturbation 
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theory is that it can also be extended to sub-horizon scales, whereas the SN 
formalism is only valid on super-horizon scales, and in some cases it has been 
shown to be numerically more efficient [107]. 

Typically non-linear effects are going to be small, given that scalar metric 
perturbations are only of order 10~^. Primordial tensor modes are smaller, 
and vector modes arc effectively zero during scalar field driven inflation. But 
the additional information available from higher-order correlations and the 
use of optimised filters for specific forms of non-Gaussianity means that there 
are already impressive constraints on the degree of primordial Gaussianity. 

Qualitatively new effects appear beyond linear order. The non-linearity of the 
field equations inevitably leads to mixing between scalar, vector and tensor 
modes and the existence of primordial density perturbations then inevitably 
generates vector and tensor modes [117,126,5,18,90]. As shown in Section 6, if 
we continue studies of scalar perturbations to higher order then the distinction 
between the different types of perturbations becomes gauge dependent and 
consistent computation will require careful (gauge invariant) definition of the 
variables being used. 
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A Definitions and notation 



A.l Notation 



The sign convention is (+++) in the classification of Misner, Thorne, and 
Wheeler [115]. 



Tensor indices: 

Greek indices, such as a, /3, . . . , /x, i/, . . . , run from to 3. Latin indices, such 
ets a,b, . . . ,i, j, . . . , run from 1 to 3, that is only over spatial dimensions. 



Spatial three- vectors are written in boldface, i.e. v = w*, whenever convenient. 
Throughout this work we use the units c — h = l. 
A. 2 Definitions 

The connection coefficient is defined as 

= 2^"^ i9ap,^i + ga^l.,p - gpix,a) ■ (A.l) 

The Riemann tensor is defined as 

pa pa pa i pa pA pa pA / \ n\ 

The Ricci tensor is a contraction of the Riemann tensor and given by 

Rfj,u — R'^ij^au I (A-3) 

and the Ricci scalar is given by contracting the Ricci tensor 

R = R^,. (A.4) 

The Einstein tensor is defined as 

G^u — R^u — ijgHvR ■ (A-5) 
The covariant derivatives are denoted by 

= V n Covariant differentiation with respect to g^v . 
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Partial derivatives are denoted by 

Symmetrisation and anti-symmetrisation are, as usual, defined as 



A. 3 Lie derivatives 

The Lie derivatives with respect to a vector field of a scalar (p, a covariant 
vector Vf^, and a covariant tensor t^i, are given by (see e.g. Ref. [160]) 

A^ = eV,A, (A.8) 

£^V^^V^,aC + Ve.^':',, (A.9) 

£^t^u = i^.,Ae^ + + ■ (A.IO) 
A. 4 Covariant derivatives 



The covariant derivatives of a scalar ip, a covariant vector V^, and a covariant 
tensor tui, are given by (e.g. [160]) 



1/ —V —r°'v 



B Components of connection coefficients and tensors 



In the following no gauge is specified, and we leave quantities undecomposed 
whenever convenient (the decomposition rules are given in Eqs. (2.4), (2.11), 
and (2.12)). 



B.l Connection coefficients 



The connection coefficients up to and including second order are 
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^ 00 



pO . 



pi 
^ 00 



IJ 



pi 



a 
a 



.a 
a 



o'^ rtik ID o/^ik jD' nr^iki 
CI 



-2Cf C(^fc - Bi I — + 01J ) + Cf (Sij- fe - Bikj) , 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



a a 
2— 

a a 



Sij + C'nj + -C2ij + 2—Cij H C2jj 
^ a Qj 



2 i^^J'i + ~ 4 (-^2j,i + -B2ij) + -Bi {Cijk,i + Ciifej - 



+01 



B 



a 



(B.5) 



"^jk, 



2 (^2fe,j + ^2i,fc ^2jA:, 



^2^1 {B'^k,3 + Bij^k) + -Bi* 



2—(f)i5kj - C[^k ~ 2— Cijfc 



+ ^.—SkjBiiCi 
a 

(B.6) 



including scalar, vector, and tensor perturbations. 



B.2 Energy-momentum tensor for N scalar fields 



The energy-momentum tensor for N scalar fields with potential U{(pi) is then 
split into background, first, and second order perturbations, using Eq. (2.1), 
as 

and we get for the components, from Eq. (5.10), at zeroth order 



^(0)0 = - (E + Uo^ , T^), = ^U" - Uo^ ,(B. 



at first order 
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a ^ 

" K 
. K 



(B.9) 



" 7 ' 



and at second order in the perturbations 



''-^(2)0 



a ^ L 



B.IO) 



K 



2 
a 



B. 3 Energy-momentum tensor for fluid 



We get for the components of the stress energy tensor, indices down, in the 
background 



Too — ~o^po ) — , Tij — a^Po , 



at first order. 



(B.ll) 



«(5Too = a2 ((5pi + 2po0i) , 

(i)(5To, = [- (po + Po) vu + Pq^h] , 



(B.12) 
(B.13) 
(B.14) 



and at second order 
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^^^SToo = [5p2 + 2po(/>2 + HiSpi + 2 (po + ^o) vikVi'' 



(B.15) 



Oi — — 



(po + -Po) (^"21 + 20it;ii + 4:CukVi 
+2vu (Spi + 5Pi) + AnukVi'' 



PoB2i + 26piB 



(B.16) 



dP2 5ij + 2PoC2ij + ^6PiCuj + (po + ^o) (vu + ^i^) + 5y) 



(B.17) 



C Geometry of spatial hypersurfaces 



C. 1 Components at first and second order of shear, expansion, and acceler- 
ation 



The calculation of the shear, defined above in Eq. (3.9), simplifies in case of 
the unit normal vector field n^, that is for = 0, 



= -nor°. - -9 Qij , 



(C.l) 



which gives (including vectors and tensors) at first order 



(C.2) 
(C.3) 



and at second order 



(j(2)(joi = 2a 



C2ij - B2{i,j) + 2BI {Cikij + Cikj,i - Cuj,k) + 201 (-Bi(i,j) - C(y) 



(C.4) 
(C.5) 



+4Cf {C[^.i - Bik,i) - 2B{ {2C^ij} - C^^ t 



(C.6) 
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The expansion is given from Eq. (3.8) in the background as 

_3a' 



(C.7) 



at first order 



and at second order 



a 



.a' 
a 



-3-01 + C,,'' - 



(C.8) 
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(Jj 



-?>-B,,Bl - ACTC[,, + ACTBu,k + 4SlC,,,f - 2BtC[ 



The acceleration is given from Eq. (3.11) at first order as 

0(1)0 = , a(i)j = 01 j , 

and at second order as 



(C.9) 



(C.IO) 



a(2)o = 2Sf0i^fe, a(2)i= 02,, + (SifcSf - 201^) . 



(C.ll) 



C.2 Curvature of spatial three-hypersurfaces 



The intrinsic curvature of spatial three-hypersurfaces is given at first and 
second order, respectively, by 



(C.12) 



-im n 
m. 



"I" ^1 y^lmn, k ' ^1 k,mn ^lmk,n ^lkn,m J 

_lO f/^^ /^j"- I n 3 k /^k /^mn \ ( (~t iq\ 



where we used 



(C.14) 
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D Governing equations 



It is often convenient to have all relevant equations available "at a single 
glance" . We therefore reproduce all governing equations given in previous sec- 
tions in this appendix together. No gauge is specified, i.e. without choosing 
a particular hypersurface or gauge restrictions, and we leave quantities unde- 
composed whenever possible (the decomposition rules are given in Eqs. (2.4), 
(2.11), and (2.12)). 



D. 1 Background 

Energy conservation for the a-fluid in the background is given from Eq. (4.21) 
as 

POa = -37^ (poa + PoJ + aQoa . (D.l) 

and the total energy conservation is then given by summing over the individual 
fluids and using Eq. (4.22) as 

p(, = -37i(po + Po) . (D.2) 

The Friedmann constraint is given from the — component Eq. (8.1) as 



The trace gives 



n" 

%-2-^87rGa^Po. (D.4) 



D.2 First order 



In this subsection we give the governing equations on large scales in the gen- 
eral case without any gauge restrictions, i.e. without choosing a particular 
hypersurface. 



D.2.1 Field equations 

The — Einstein equation is given from Eq. (8.1) as 

m {Hh + ^0 - V' (^1 + Hai) = -AnGa'^Spi . (D.5) 
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The — i Einstein equation is 



(D.6) 



the off-trace is 



and the trace is 



(D.7) 



jjj'l + 2ni^[ + 7^0; + (^2^ - ^ j 0, = AnGa^ [sP^ + ^V'n) . (D.8) 



D.2.2 Energy-momentum conservation 

Energy and momentum conservation of the ct-fiuid is given from Eq. (4.21) at 
first order as 

V„+37i(5p„+5P«)-3 (p„ + P„) V;'+a-i(p„+P„)V' {V^ + a) = aQ^cf>+adQa , 

(D.9) 

The momentum conservation equation of the a-fiuid is 



aQc 



pa + Pa 



[1 + cl) - mcl 



6P^ + -v^n, - Q^V - 



Pa + Pa 

(D.IO) 

Total energy and momentum conservation follows from the above, by summing 
over all individual fluids and using Eqs. (4.22), and is given by 



= 0. 



Vi + 37i ((5pi + (5Pi) + (po + Po) ((7i + ^;l + Pi) - 3V^;] = , (D.ll) 

[(po + Po) (^1 + Si)]' + (po + Po) m (^1 + Si) + 0i] + ^Pi + ^ v'n = . 

(D.12) 



D.3 Second order 



D.3.1 Energy-momentum conservation 

In the multi-fluid case, energy conservation of the a-fluid is given from Eq. (4.21) 
at second order as 
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+ 2 {poa + Poa) [{v^a)k + ^Ife) {^V^a)k + ^Ifc) + 4HVi%) (i;i(a)fe + ^lik) 
-4 (poa + ^la) Ci.^.Cj^' + 2 (poa + ^la) (S^^L^l,^ + (pl^^Via) 

+27ri(„),, (Cf^- - 2ncr) + 27:1l^v^^^),,l + 2^;f,7r,(J),,, 
+2 (po« + Poa) (-3^1,; + V'^i,;) 

= a|(5g2a + 20i5Qi« + Qoa (</'2 - 01^ + -yife^^i) + ^/i(a)fe^^i I ■ (D-13) 



For a single fluid we flnd by summing over the individual fluids and using 
Eq. (4.22) 



V2 + 37i (5p2 + 5P2) + (po + Po) (-3^/'2 + V'£;^ + V^V2) 

+2 {6pi + 5Pi) (-3^; + V^E[ + + 2 {6pi + 5Pi)_, + 2 (p[, + P^) t;^' + P^) 

+2 (po + Po) + B[„) {2v,k + Pi^) + mv', (vik + ^i^)] - 4 (po + Po) C[,fl' 

+2 (po + Po) {2v^<Pi,k + ^i^^v,) + 27ri,,- (cf^' - 27iCi^) + 27r^,W,i + ^v^tt^J^ 

+2 (po + Po) {-3^1,1 + V^E,,i) v[^0. (D.14) 
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List of Symbols 



a Scale factor 

Acceleration 
Cg Adiabatic sound speed 
Cq, Adiabatic sound speed of the a component 
/(J) Momentum transfer perturbation of Jth component 
Qfj^i, Metric tensor 
hij Tensor metric perturbation 
V,j,u Projection tensor, V^,, = g,j,i, + n^n„ 

Comoving wave-vector 
k Comoving wavenumber k'^ = k^ki 
Mpi Planck mass = 
mpi Reduced Planck mass = (SttG)"^ 

Unit time-like vector field 
Ug Spectral index of curvature perturbations 
Anx Scale dependence of perturbation spectrum of a quantity X 

momentum current density (in rest frame) 

Wave vector 
ds Infinitesimal line element 
t Coordinate time 
iiiJ- 4-velocity 

V Scalar velocity perturbation 

Va Scalar velocity perturbation of a-component 

vj Rescaled perturbation in Jth orthogonal field, a5aj 

Vector velocity perturbation 
X* Spatial coordinate 
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B Shift (scalar metric perturbation) 

E Scalar spatial metric perturbation 

Fi Vector spatial metric perturbation 

G Newton's constant 

G^i, Einstein tensor 

H Hubble parameter, H = ^ 

H^P Hankel function of the first kind of degree u 

K^i, Extrinsic curvature 

L Lagrangian 

jC Lagrangian density 

N Number of e-folds (integrated expansion) 

Vfj,u Projection tensor, V^i, = g^^ + n^n^, 

P Pressure 

Pa Pressure of the a-component 

Energy transfer parameter of the Jth component 

Q^j) Energy momentum four vector of the Jth component 

Vx Power spectrum of a quantity X 

R Ricci scalar 

Rfj,v Ricci tensor 

^^^R Intrinsic spatial curvature of three-hypersurface 

TZ Curvature perturbation in comoving gauge 

S Action 

Si Vector metric perturbation 

Sij Entropy perturbation 

T^i, Energy momentum tensor 

U {(f) Potential of scalar field 
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a Arbitrary scalar function (temporal gauge function) 

/3 Arbitrary scalar function (spatial gauge function) 

Kronecker delta 

77 Conformal time, adr] = dt 

C Curvature perturbation on uniform density hypersurface 

Y Arbitrary divergence-free vector function (spatial gauge function) 

7ife Metric tensor on spatial 3- hypersurface 

K Curvature of background spacetime 

A comoving wave length k — ^ 

'^phys physical wave length Aphys = aX 

cUfj^i, Vorticity tensor 

TT* Anisotropic stress vector 

TT^j^ Anisotropic stress tensor 

tensor^i Teusorial anisotropic stress tensor 

if Scalar field 

(fii one of multiple scalar fields 

(f) Lapse function (scalar metric perturbation) 

■0 Curvature perturbation (scalar metric perturbation) 

p Energy density 

Pa Energy density of a-component 

a Shear scalar 

afj,i, Shear tensor 

T proper time, o?t^ = ds"^ 

9 Expansion 

Arbitrary vector valued function 

r(a;) Gamma function 

n Scalar anisotropic stress tensor 

$ Bardeen potential (lapse function in longitudinal gauge) 

^ Bardeen potential (curvature perturbation in longitudinal gauge) 



85 



Dimensions 



It can be useful as a quick check of the vahdity of an equation or expression, 
particularly for the large expressions at second order, to check that all terms 
have the correct dimensions. 



Quantity 


Dimension 


[a] 


1 


M 


rp-l 


[H] 




[U] 


rp—4: 




rp-2 


[G] 


rp2 


[v] 


T 




T 




T 


[0] 


1 




1 


[B] 


T 


[E] 






1 


[7^] 


1 


[C] 


1 


[p] 


rp—i 



Note, in geometric units T = L. 
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